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ION EXTRACTION CAPABILITIES OF 
TWO-GRID ACCELERATOR SYSTEMS 

D. C, Rovang 

An experimental investigation of the ion extraction capabilities of 
two-grid accelerator systems was conducted during this grant period. This 
work. is described in Ref. I so it will not be reproduced here._Jhe body of 
experimental data presented in Ref. 1 facilitates the selection of the 
accelerator system geometries and operating parameters necessary to maximize 
extracted ion current. Reductions in screen hole diameter to 0.5 mm are 
shown to induce no dramatic changes in impingement-limited perveance. Ra- 
sul ts_obtained at small grid separation ratios suggest a new grid operating 
condition where high beam current per hale levels are achieved atJL specified 
net accelerating voltage. This operating condition is realized at a relatively 
"*■ high, optimum ratio of net-to-total accelerating voltage ratio. Electron 
backstreaming and electrical breakdown characteristics of two-grid systems 
are also given. 

ELECTROTHERMAL RAMJET RESEARCH 
B. D. Shaw 

Extensive theoretical analysis of the electrothermal ramjet^ concept 
was completed during this grant period. This work describes the effects of 
ionization and dissociation of the propellant on the performance of the ramjet. 
The effects of ionization are shown to be small so they can be neglected. 
Dissociation phenomena ace modelled using equilibrium, frozen flow and 
finite reaction rate approaches and the equilibrium approach is shown to 



SIMPLE PERFORMANCE-MODEL FOR RING AND LINE CUSP ION THRUSTERS ..- 

John R. Brophy 

Introduction 

In order to .improve thruster performance it is desirable to have a 
theoretical model which describes the effects of thruster design variables 
and operating parameters on performance. Theoretical models have been 
developed in the past to predict the performance of both low magnetic field 
strength (Kaufman type) thrusters 3 and high field strength cusped thrusters. 4 
However, these models typically consist of a complex set of equations which 
must be solved iteratively by a computer. This complexity makes it diffi- 
cult to identify the dominate mechanisms controlling thruster performance. 
Earlier work. (on Kaufman type thrusters) by Masek 5 and Knauer 6 described 
thruster performance physically in terms of the average energy expended to 
produce ions in -the discharge chamber plasma. This work resulted in an 
increased understanding of some of the phenomena affecting thruster per- 
formance, but. stopped short of providing the capability to predict the 
performance. 

In this paper, a model is developed for high field strength cusped 
thrusters, which results directly in a single equation for the energy cost _ 
of a beam ion (eV/beam ion) as a function of the propellant utilization. 

This is made possible by formulating the model . in terms of the average 
discharge plasma ion energy cost and the fraction of ions produced which 
are extracted into the beam. A key feature of the model is that it makes 
it possible for one to calculate readily the average energy required to 
produce an ion in the discharge chamber plasma as a function of . the pro- 
pellant mass flow rate and propellant utilization. Comparison of calculated 
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and measured values of the average plasma 1on_energy cost are made for 
variations In propellant flow rate,, propellant type, discharge voltage, 
propellant utilization, Accelerator grid open area fraction and beam diam- 
eter. The effect of the variation In each of these parameters on. the ratio 
of beam current to ion current produced (extracted ion fraction) is. also 
investigated. 

Finally, the model is exercised to indicate the effects on thruster 
performance due to variations in such parameters as; propellant flow rate, 
propellant type, discharge voltage, extracted ion fraction, accelerator 
grid transparency to neutrals and the containment of primary electrons. The 
potential for improvements in thruster design suggested by the model are 
discussed. 


Theoretical Development 


Beam Ion Energy Cost 

The specific discharge power, or energy cost per beam ion is defined as, 

e B = ^ J D " J B^ V J B* C eV / beam i or 0 (1) 

where the symbols are defined in Appendix A. The beam current (Jg) in Eq. 1 
is subtracted from the discharge current ( Jq) so that the energy that goes into 
accelerating the beam ions through the discharge voltage is not charged to 
the energy cost per bearn ion. 

In a similar manner -the average energy expended in creating ions in 
the discharge chamber plasma may be defined as, 

e p B C J o “ ( J C + V J p ’ CeV/plasma ion]. 


( 2 ) 
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By analogy to_Eq. 1,. the + Jg" term Is. subtracted from the discharge 
current so that , the energy that goes Into accelerating these 1ons_out of 
the discharge chamber plasma Into the-chaniber walls or the beam Is not 
Included In the Ion production cost. Rearranging Eq. 2 yields, 


e 


P 


(J D~ J B ) 



(3) 


Defining the .fractions of ion current produced that go into the beam and to 
cathode potential surfaces as. 


f B * V J p and f C * J C /J p ’ < 4 > 

respectively, and using Eq. 1 in Eq. 3 yields, 

e p = £ B f B “'Vd * 

Solving this equation for gives the result, 

e B = e p^ f B + ' f C V D //f B * 

This. equation describes the beam ion energy cost as a function of the plasma 
ion energy cost (e )., the extracted ion fraction (fg), the fraction of ion 
current to cathode potential surfaces (f^) and the discharge voltage (V D ). 

The first term on the right-hand side of Eq. 6 represents the energy 
loss associated with producing ions in the discharge chamber and extracting 
only a fraction of them into the. beam. Ions which are not extracted into 
the beam go to the walls of the discharge chamber where they recombine. The 
resulting atoms must then be re-ionized before they can contribute to the 
beam current. Thus the. factor 1/fg may be interpreted as the average number 
of times that a beam ion undergoes ionization before being extracted into 
the beam. 
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The second tern on. the right-hand side of Eq. 6 represents the energy 
wasted In accelerating plasma Ions Into Interior cathode potent1al_surfaces. 
TMs process results In both a discharge energy loss and In the sputter 
erosion of these surfaces. 

To generate performance curves using Eq. 6 one must be able to specify 
the behavior of each of the terms on the right-hand side of Eq. 6 as a 
function of the propellant utilization. 

Plasma Ion Energy Cost 

The plasma ion energy cost parameter, 6p , appearing in Eq. 6 and 
defined by Eq. 2 reflects all mechanisms of energy loss from the discharge 
chamber except for the acceleration of ions out of the plasma. through the 
discharge voltage. Specifically Ep includes energy losses due to the follow- 
ing mechanisms; direct primary electron loss to the anode, Maxwellian 
electron collection by the anode, excitation of neutral atoms, excitation, 
of ionic states (which will be neglected) and hollow cathode operation. To 
derive an expression for the plasma ion energy cost as a function .of the. 
propellant utilization an energy balance is made on the discharge chamber 
plasma represented in Fig,, l . The primary electrons are assumed to be 
accelerated from a cathode plasma potential that is V Q volts above cathode, 
potential to the potential of the bulk plasma which is assumed to be near 
that of the anode (i.e. .at the discharge voltage (V Q ) above cathode potential). 
Under this assumption as well as the assumption that only the discharge 
power supply is used to sustain the discharge, .energy is supplied to the 
plasma by the primary electrons at the. rate,. J £ (V D -V C ). The "missing- 
energy J E V c .1s used to operate the hollow cathode. Energy is lost from the 
plasma_primarily by the flux of four types of energy carriers across the 
plasma boundaries; ions, photons (emitted by de-excitation of excited pro- 
pellant atoms), Maxwellian electrons, and primary electrons. The ions and 
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photons. are lost to. all interior thruster surfaces whereas the Maxwellian 
and primary electrons are assumed to be lost to the anode surfaces, only. 

In steady state the rate Of. energy supplied to the plasma must be equal to 
the rate at which it is lost, thus, 


V V D -V = V + + £ J J U J + J M e M + J L< V D - V C> • (7) 

J 

Dividing Eq. 7 by the ion production current (J ) and recognizing that the 
emission current (0^) is related to the discharge current by, 

J E * J D - tJ C + J B> • W 


allows Eq. 7 to be written as, 

l 

i 




■ V 


, J M e M ^ J L V D °E V C 

"V "V "5" 


J L V C 


(9) 


4* h 

where Eq, 2 has been used. The rate at which the j excited state is ... 

produced is given by 


Jj ■ en 0 n e <0 J V *P 


(JO) 


x U 

where <OjV e > represents the product of j excitation collision cross section 
and the electron velocity averaged over the electron speed distribution. 
Similarly, the ion production current is given by, 


J p = e Ve <0 + v e > *p ’ 


( 11 ) 


Substituting Eqs. 10 and 11 into 9 yields, 


. .« + *£s + Vd t Vc . 

p o J n TT J n 

P P P P 


( 12 ) 


where , 
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e o a u + + j 




> u. 
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The third term on the right hand side of Eq. 12 may be written. 



(14) 


where the last step was made using Eqs. .2 and 8. The ratio J L /J E is simply 
the fraction of primary electrons emitted by the cathode which are collected 
by the anode before having any inelastic collisions. This fraction may be 
given by the survival equation 7 as, 


\ _ 'Vo ‘e 

J E 


(15) 


where is the average distance a primary electron would travel in the dis- 
charge chamber before being collected by the anode - assuming it had no inelastic ... 
collisions and o Q isthe total inelastic collision cross section. Combining Eqs. 12, 
14 and 15 yields, 


. e + Vm . 'We . V C 

e ° V p + H 


\ "°o n o £ e 

e P v D e 


(16) 


The current of Maxwellian electrons to. the anode may be given as the sum of 
the secondary electrons liberated in the ionization process and the therm- 
al ized primary electrons, thus, 


J M “ J p + < J E ' J L } = J p + J E {1 ‘ °L /J E ) 


(17) 


Using Eqs. 15 and 17 in 16 yields 


( 


e n = G n + J n + Jr(l ’ e 


'o o*e 


> hr 


+ e. 


-a.h A i 

* 0 0 6 , /. 
e + (1 - e 


-^oVex V C 

' V~ 
V D 


) 


( 18 ) 
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Finally, using Eqs. 2 and 8 and solving Eq. 18 for e results in, 

r 


e o + fc M 


1 - 


V " '+ "e ' 

V C M 


1 - e 


■o 

o o e 


-1 


( 19 ) 


The neutral density, n Q , may be. expressed in terms of the propellant flow 
rate and propellant utilization by equating the rate at which propellant 
enters and leaves the discharge chamber,, i.e., 


m - J 3 + n Q (20) 

• • 

where m and n Q are in units of equivalent amperes. The neutral flow rate 
from the thruster may be given as 


"o " T "o ev o A g *0 • 


Combining Eqs. 20 and 21 yields, 

_ 4 m (l-n u ) 

n o " ev_ A„. 

o g. y o 

where n u = Og/m was used. Thus Eq . 19 may be written, 

-t 


* 

C * E 

P P 


1 - e 


•CAl-nJ 


where, 


4 a JL 

r - o e 

o ■ ev ' A' r 6 * 

o q v o 


and, 


G * E o *; e M 

- — o — - 

V D 


( 21 ) 


( 22 ) 


(23) . 

(24) 


1 


(25) 
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Equation 23 provides a very simple method, for calculating the plasma ion- 
energy cost as a function of the propellant utilization. Experimental re- 
sults, which will be presented in the next section^ indicate that under many 
conditions the parameters C Q and may be taken to be independent of the 
propellant utilization. Substitution of Eq. 23 into Eq. 6 yields a single 
equation describing the performance of a given thruster design. For design 
purposes the parameters fg and f^ in addition to C Q and e* may be ..taken to 
be independent of the utilization and flow rate. These parameters do, how- 
ever, depend strongly on the thruster design. Indeed, these. four. parameters 
determine the performance of a given thruster design. 

The parameter C Q depends on the quality of the primary electron contain 
ment (through & e ) , the quality of the containment of neutrals (through Ag, 

4*0 and V’ - and the Propel! ant.. type (through a Q and v Q ). Recall that the 
primary electron containment length may be interperted as the average 
distance a primary electron would travel in the discharge chamber before 
being lost to the anode - assuming.it had no inelastic collisions. Magnetic 
fields in all discharge chamber designs serve the function of increasing 
this length. Although an effective means of determining a remains to be 
developed it is believed that this parameter is a function primarily of the 
thruster geometry, magnetic field configuration and cathode location. 

Through the parameter C Q , Eq. 23 suggests that the plasma ion energy 

cost should depend on: the propellant; since different propel 1 ants. will 

have different values of the total inelastic collision cross section (a ) 

v o' 

and different atomic masses, which implies different neutral velocities; the 
wall temperature, giving rise to different neutral velocities (v Q ); the 
transparency of the grids to neutral atoms U Q ); the area through, which. the 
beam is extracted (Ag) and finally, the discharge voltage which determines 
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the primary electron energy and thus affects the value of the cross .section 

V 

The parameter e p depends on a number of energy loss mechanisms including 
the relative amount of energy expended in excitations compared to ionizations 
of neutral atoms through e Q , the. average energy of the Maxwellian electrons 
wmch leave the plasma, e m , and the efficiency with which the hollow cathode 
operates. The cathode efficiency is reflected in the value of V r which 
represents plasma potential from which the electrons are supplied. Inefficient 
cathode operation results in high values of and corresponding poor overall 
thrus-.er performance. For thermionic cathodes = 0. , however, additional, 
heater power must be supplied to effect their .operation. The parameter.^, 
defined by Eq. 13, in general depends on the electron energy distribution 
and may be calculated in the manner discussed by Dugan and Sovie? In spite of 
the complexity associated with the computation of the parameter e* it will 
be shown later that for design purposes it may be taken to be constant for a 
given discharge chamber design, propellant type and discharge voltage.. Thus 
the model does not necessarily require a detailed knowledge of the electron 
energy distribution as a function of the propellant utilization to evaluate 
different thruster designs.. 

Although the model cannot yet be used to predict the performance of 
completely new thruster designs it provides a clear physical picture of the 
phenomena affecting the performance.. In order .to test several key features 

of the model the experimental investigation described in subsequent sections 
was undertaken. 
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Experimental Apparatus 

For this investigation the ion source shown schematically in Fig. 2 
was designed and built. This source normally produces a 12 cm dia. ion .beam 
and provides the capability for measuring the distribution of ion currents 

to the beam, screen grid and internal, thruster surfaces (with the exception 
of the anode). 

The magnetic field for this experimental ion source is established 
through the use of an electromagnet located on the upstream centerline of 
the discharge chamber and a number of 1.1 cm x 1.3 cm x 0.5 cm samarium 
cobalt permanent magnets.. These permanent magnets are arranged end-to-end 
to form ring magnets of alternate polarity in the manner suggested by Fig. 2. 

I he flux density at the surface of the magnets is 0.27T and the magnets are 
attached to the steel discharge. chamber housing by their own magnetic attrac- 
tion. This arrangement allows the ion source magnetic field configuration 
to be altered quickly and easily by simply adding, removing or changing the 
position of the magnets. Although many different configurations were tested 
the results obtained were all similar, thus only those obtained using the 
configurations shown in Figs. 2 and 3 will be presented. For the configuration 
shown in Fig. 2 the upstream magnet ring is covered with a strip of 0.13 nrn 
thick steel insulated from the magnets themselves by a strip of 0.25 mm thick 
mica. This is done so that the surface of this strip can be maintained at anode 
potential while the rest of the thruster body is at or negative.of cathode poten- 
tial. The downstream magnet ring is uncovered. The magnetic flux density at the 
surface of the electromagnet can be adjusted from zero to approximately 0.2 T 
by adjusting the magnet current from zero to 124 A. 

The main discharge chamber cathode consists of seven 0.25 cm dia. 
tungsten wires connected in parallel and supported by two support posts that 
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are electrically .Isolated from, the thruster body. Each cathode wire Is 
approximately 2.8 cm long .so. the total, cathode length exposed to the plasma 
Is about 19.6 cm. These seven short wires ,.in parallel are used to minimize 
the voltage drop across the cathode, A voltage drop less than 3 v at the 
maximum heater current was achieved with this system. It is noted that a 
small voltage drop across the cathode results in a primary electron energy 
distribution that more closely resembles the monoenergetic distribution 
produced by a hollow cathode. The cathode wires were heated using direct 
currents in the range 6 to 8 A per wire. Tests were conducted using both 
argon and krypton propellants. Discharge voltages were varied from 30 to 
50 v for argon and 20 to 40 v for krypton. The discharge current was adjusted 
through the range of 0.5 to 5 A for both propellants by controlling the heater 
current through the refractory cathode wires. 

Two ion accelerator systems were used in this study. _The first accel- - 

erator system consisted of a set of dished small hole accelerator grids 

(SHAG) with a cold grid separation of 0.75 mm and screen and accelerator grid 

physical open area fractions of 0.68 and 0.3CL, respectively. The second 

system consisted of a set of dished large hole accelerator grids (LHAG) with 

a cold grid separation of 0.75 mm, and screen and accelerator grid physical 

open area fractions of 0.68 and 0.57, respectively. Both accelerator systems 

were normally masked to produce a 12 cm diameter beam. One series of tests 

was conducted, however, with the SHAG set masked to produce a 6 cm diameter 

beam. For the 12 cm dla, beam tests flow rates for both argon and krypton 

were varied from. 500. to 1500 mAeq. For the 6 cm dia, beam test th? flow 

rates were varied from 125 to 500 mAe.q, All tests were conducted in a 1.2 m 

-6 

dia. x 4.6 jn long vacuum test facility. Tank pressures ranged from *2 x 10 
Torr with no flow to ^3 x 10"^ Torr at a flow rate of 1500 mAeq of krypton. 
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A second series of tests were performed on the line. cusp thruster con- 
figuration shown schematically In Fig. 3. The magnetic field for the 
configuration is established using eight, rows of samarium cobalt permanent 
magnets of alternating polarity positioned along the. cylindrical discharge 
chamber side wall. Six of the magnetic rows were maintained at anode po- 
tential while the other two (diametrically opposite) rows were biased 
negative of cathode potential along with the back surface, side wall and 
screen grid. -Tests w ere conducted over a range of argon propellant flow 
rates of. 600 to 1600 mAeq, discharge voltages of 40 to 50 v and discharge 
currents of 0.5 to 4.0 A. 


Procedure 

The following set of experiments was designed to test the ability of 
the model just developed to predict the functional dependence of the 
parameter £p. The model predicts that the plasma ion energy cost should 
behave according to Eq. 23 with C Q and Ep given by Eqs. 24 and 25, 
respectively. The value of ^ may be determined experimentally through the 
use of Eq. 2 provided the total ion production current (Jp) and the ion 
currents, to the beam (Jg) and to cathode potential surfaces (J r ) are 
measured. The total ion production current may be given as the sum of the 
ion currents leaving the plasma, i.e., 

J P " J 6 + J C + J A • (26) 

For these tests the thruster configuration of Fig. 2 was operated with 
only the upstream magnetic ring at anode potential. All other interior 
discharge chamber surfaces (with the exception of. the cathode support 
posts) were biased 30 v negative of cathode potential to repel the 
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discharge chamber. electrons and facilitate the measurement of the Ion cur- 
rent to those surfaces. Since the only anode surface occurs at a magnetic 
field cusp the effective area for Ion loss to this surface Is expected to 
be less than the physical area. 9 ' 12 Rough calculations Indicate that the 
ion current to the anode, with this configuration should be less than a few 
percent of the total production current. Thus J may_.be approximated 

r 

as the sum of the -ion currents to the beam (including the impingement cur- 
rent) and to the negatively biased discharge chamber surfaces (including 
the screen grid). 

Data consisting of the beam current, propellant flow rate, propellant 
utilization and total ion production current were collected over the range 
of operating conditions discussed earlier with the electromagnetic current 
held constant at 57 A for all cases. Thus the plasma ion energy cost was 
determined from Eq. 2 for a wide range of operating conditions, in addi- 
tion, the extracted ion fraction, f g , given by Eq. 4, was also measured 
over this same range of conditions. 

Results and Discussion 

Plasma Ion Energy Cost 

Measurement of the plasma ion energy cost for operation of the ring cusp 

thruster with argon propellant at a 50 v discharge voltage over the range of 

neutral flow rates from 500 to 1500 mAeq yielded the results shown in Fig. 4. 
Here the measured values of e p are plotted as a function of the quantity 
m(l - n u ) as suggested by Eq. 23. Since this quantity. is related to the 
neutral density in the discharge chamber [Eq. 22] it will be referred to as 
the neutral density parameter. The solid line in Fig. 4 is the curve given 

by Eq. 23, where the parameters C Q and s* have been selected to give the best 
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fit to the data. The parameter t* was taken to be the value of c measured 
at-large values of the neutral density parameter. The reason for this can 
be understood by considering Eq. 23, which shows that when C Q m(l - n u ) Is. 
large then the exponential tern, Is small compared to unity and one obtains 
o p * £ p . Having established the value of c* the value of C 0 is varied until 
the best fit Is obtained. The agreement between the functional form of 
Eq. 23 and the experimental dataJs seen to bex,uite good. This Indicates 

that the parameters C 0 and E * may indeed be taken to be independent of the 
neutral density parameter. 

Ue have now established a value of c„ . 3.1 Aeq' 1 . which is applicable 

to the ion source of Fig. 2 operating at the conditions defined In the legend 

for Fig. 4. New values of C 0 applicable to other operating conditions may 

now.be calculated from this value using Eq. 24. For example, changing grid 

sets from SHAG to LHAG should change the value of C Q through the P ara re ter 

4 0 which is the transparency of the grids to neutral atoms. This effective 

transparency parameter may be calculated for each grid set according to the 
equation, 

+ s <f» a 

*° = (27) 

where ♦, and ^ are the modified transparencies for the screen and accelera- 
tor grid, respecti vely. These modified transparencies may be calculated as 
the physical open area fraction of a grid times the appropriate clausing 
factor. For the two grid sets used in this study we have, 




nr = 0J6 » and * 0 “ 0.27 

SHAG °LHAG 


( 28 ) 


Thus the new value of C Q applicable. to the LHAG optics with all other 

i I* . * . . . - . 


conditions held constant is given by, 
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which yields (C q ) lh ^q =1.8 Aeq 


-1 


The measured values of e obtained under the same set of .conditions 

r 

defined in. the legend of Fig. 4 except for the change in optics from SHAG 

to LHAG yielded the results shown by the data points in Fig. 5. The solid 

line is the prediction of ..the model based on the value of C Q calculated from 

★ ★ 

Eq. 29. The value of e was held constant at e = 57 eV since changing the 

r K 

optics should not affect this parameter. Clearly the model correctly predicts 
the variation in the plasma ion energy cost. . 

The same procedure of calculating a new value of C Q from the old value 
according to Eq. 24 was followed for the analysis of the data displayed in 
Figs. 6 through 9. For. the data in Fig. 6 the thruster was operated with 
krypton propellant and SHAG optics at a discharge voltage of 40 v. 

The corresponding value. of C Q was calculated using the value of C 0 
obtained from Fig. 4 and the equation, 


(CJ 


Kr 



(CJ 


Ar 


(30) 


where. both the change in propellant properties and the change in discharge 
voltage must be accounted for. That is, (o Q ) Kr 1s the tota ^ inelastic 
collision cross section for 40 eV primary electron - krypton atom collisions, 
whereas (o Q ) Ar refers, in this case, to 50 eV primary electron - argon atom 
collisions. The cross section data needed in this equation were obtained 


from de Heer, et. al . 
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The new value of C Q calculated from Eq. 30 was _ 


C q = 5.7 Aeq ...... and the corresponding prediction of the model is compared 
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to the measured values In Fig. 6. A new value of c* was also used for the 

P 

model In Fig, 6 since this parameter. should be a function of both the dis- 
charge voltage and the. propellant type. Again this new e* was obtained 
according to the procedure outlined earlier. 

Continuing the same procedure, the prediction of the model based on a 
calculated new value of C Q is compared to the experimental results for 
operation with LHAG optics, krypton propellant and V Q = 40 v, in Fig. 7. 

Again, the agreement between the predicted curve and the experimental data 
is seen to be good. .... 

From Eq. 24 it is seen that C Q depends inversely on the area of the 
grids through which the beam is extracted, A g . This parameter may be varied 
without changing the discharge chamber diameter by masking down the screen- 
grid to produce ion beams of different cross sectional areas. In this case 
the screen grid for the SHAG optics set was masked down from a beam diameter 
of 12 cm to one of 6 cm.. This fourfold reduction in beam area should pro- 
duce a corresponding fourfold increase in the parameter C Q . Therefore, 
taking C Q from Fig, 6 and multiplying by four yields the new value of 
C Q = 22.8 Aeq 1 . The model prediction using this value of C Q is compared 
to the measured values of Sp in Fig. 8. Remarkably, the agreement between 

the model and the experiment is excellent. Similar agreement was obtained 

for operation with argon using the masked down grid set. 

Next, the dependency of the plasma ion energy cost on discharge voltage 
for a given propellant is investigated. For operation with argon at V D = 40 v 
a new value of C Q was calculated from the value given in Fig. .4 which was 
obtained at Vq - 50 v. Again the model agrees well with the measured values 
as shown in Fig. 9, provided a new value of e* = 60 eV is used. 

For operation at V Q =30 v with argon, however, the situation is quite 
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different as seen in Fig, 10. Here a systematic difference in the data taken 
at different flow rates is observed. Clearly a single equation such as 
Eq. 24 is not sufficient to explain this behavior if C Q and e* are taken to 
be independent of the neutral flow rate and the utilization. Since the 
neutral density parameter may be expressed as m - J B we see that a. given 
value of this parameter obtained at different flow rates corresponds to 
different beam currents. Further, the different beam currents may be roughly 
translated into different plasma densities. Thus a possible theoretical 
interpretation of the data in Fig. 10 might be found by looking for C Q and/or 
£ p t0 be Unctions of the plasma density at low discharge voltages. 

For krypton propellant, this separation with flow rate was found to occur 
at a discharge voltage (V Q ) of 20 v rather than at the 30 v value observed 
for argon. This indicates that the separation phenomena depends on the 

properties of the iens (or neutrals) in the discharge chamber as well as 
the discharge voltage. 

Experiments performed on the line cusp configuration of Fig. 3 indicate 
that the model is applicable to line cusp as well as ring cusp thruster 
designs. An example of the results obtained with the line cusp discharge 
chamber is given in Fig. 11 for operation with argon propellant at. a dis- 
charge voltage of 50 v. Good agreement between the model and the experiment 
is obtained for values of C Q = 1.6 Aeq -1 and e* = 50 eV. 

The results discussed so far indicate that the model for the plasma 

ion energy cost correctly describes the behavior of. this parameter over a 

wide range of operating conditions. It also indicates that, with the ex- 
ception of operation at low discharge voltages, the parameters C and 
* ^ 
e p may be taken to be independent of the neutral flow rate and beam current. 

Physically, Eq. 24 describes the plasma ion energy cost in terms of the 
loss of primary electrons to the anode. At high values of the neutral 
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density parameter the neutral density in the discharge chamber is- large and 
the probability Is high that all the primary electrons will undergo inelastic 
collisions with, neutral .atoms and none will be lost directly to the anode. 

In this case the discharge chamber will be producing ions for the minimum 
energy cost, e p . -As the beam current is increased (for a constant flow rate) 
the propellant utilization increases causing the neutral density parameter and 
therefore the neutral density to decrease (see Eqs. 20-22). The decrease in 
neutral density increases the likelihood of a primary electron reaching the 
anode- without first losing its energy. This direct loss of primary electron _ 
energy increases the overall plasma ion energy cost according to Eq. 24 and 
consequently increases the beam ion energy cost according to Eq. 6. Indeed, 
the shape of the performance curve is largely determined by the direct loss 
of primary electrons. Thus, any design change which decreases the likelihood 
of the direct loss of primary electrons (without decreasing the extracted 
ion fraction) should improve, the thruster's performance.. The probability of 
direct primary electron loss is determined through the parameter C Q . This 
parameter may be increased by either increasing which makes it harder 
for primary electrons to escape the plasma, or by making it harder for neutral 
atoms to escape the discharge chamber. 

Now, because the plasma ion energy cost depends, on the neutral density 
parameter and hence on the neutral density itself we are in a position to 
understand the difference in discharge chamber performance observed for 
thruster operation with and without ion beam extraction. .It has been observed 
that the performance (eV/beam ion) extrapolated from thruster operation 
without beam extraction is generally significantly better than the performance 
measured with beam extraction. 16 ’ 16 This may be explained by noting that 
thruster operation without beam extraction should be characterized by higher 
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discharge chamber neutral densities than operation with beam extraction at 
the same flow rate. This must be so since without beam extraction most of 
the propellant leaves the discharge chamber as neutrals. For the case with 
ion extraction most of the. propellant leaves as ions. Consequently, fewer 
neutrals leak out through the grids with ion extraction as compared to opera- 
tion at the same flow rate without ion extraction. This v.ranslates into 
lower neutral densities according to Eq. 21, for the beam extraction case. 
Further, Eq. 23 written as 

* 

= e 

P P 

indicates that the plasma ion energy cost should be higher for smaller 
values of n Q . Thus, the performance with beam extraction should be poorer 
than that extrapolated from data obtained without beam extraction. Finally, 
it is noted that careful application of the performance model developed here 
should allow more meaningful extrapolation of data taken without beam ex-, 
traction to operation with beam extraction. 

Extracted Ion Fraction 

Aside from the plasma ion energy cost, the other parameter which has a 
major affect on. thruster performance is the extracted ion fraction (f g in 
Eq. 6). It is of little use to create ions efficiently in the discharge 
chamber if only a small fraction are extracted into the beam. The extracted 
ion fraction, which is defined as the beam current divided by the total ion 
production current, was measured over the same set of operating conditions 
described earlier. The resulting data was found not to correlate well with 
the propellant utilization, where a systematic difference between data taken 
at-jdifferent flow rates was observed for all cases. The extracted ion 
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fraction data was found, however, to. .correlate well with the neutral density 
parameter .for a reason that has not yet been justified physically... . 

Values of fg are presented in Fig. 12 as a function of the neutral 
density parameter for operation with argon propellant.at discharge voltages 
of both 30 and 50 v over a range of flow rates from 500 to 1500 mAeq, 

Several things are evident from these data. First, the degree of correlation 
for the data taken at different flow rates and a fixed discharge voltagejs 
quite good.. Second, the extracted ion fraction is greater at the lower, dis- 
charge voltage. This was true for all ..sets of data taken. The data taken 
under the same conditions but at Vg = 40 v falls between the two. curves shown 
in Fig. 12. Third, f g is seen to increase slightly with decreasing values 
of the neutral density parameter. . Finally, the data taken at Vp = 30. v 
corresponds to the data in Fig. . 10, where we see that the extracted ion 
fraction still correlates with the neutral density parameter at this low 
discharge voltage even if the plasma ion energy cost does not. 

The extracted ion fractions for argon propellant are compared to those 
for krypton for Vp = 40 v in Fig. 13. The values of fg for argon are seen 
to be generally higher than those for krypton. 

It should be pointed out that none of the effects exhibited in Figs. 12 
and 13 are understood at this time. Earlier studies indicate that fg is 

strongly dependent on the magnetic field configuration. * However, since 

f D is not a strong function of the neutral density parameter and. tends to 

D 

increase only slightly with decreasing values of this parameter, a simple 
design approach would be to take fg to be a constant for a given discharge 
chamber design (magnetic field configuration), propellant, and discharge 
voltage. 


EXTRACTED ION FRACTION, 


OR1GSNM. FW* " 
OF POOR QUAulT 


OPEN SYMBOLS 

Vq=50V 

m 

(mAeq) 

CLOSECL SYMBOLS 

V d *30V 

0 

500 



□ 

750 

ARGON 


A 

1000- 

SHAG OPTICS 


O 

1500 


NEUTRAL DENSITY PARAMETER, m(l-« ), [Aeq] 


Pin 12. Extracted Ion Fraction Results 


a 





tXTRACTED ION 


0.4 



OPEN SYMBOLS 

ARGON 


m ( mAeq) 

0.3 

CLOSED SYMBOLS 

krypton 

O 

500 




□ 

750 

0 2 

SHAG OPTICS . 


A 

1000 


> 

o 

<*• 

H 


O 

1500 


0.1 }- 


® 0.2 0.4 0.6 0.8 | # o |2 

NEUTRAL DENSITY PARAMETER, rh(l-^ u ), [ A ©q j 


Fig. 13. Extracted Ion Fraction Results 




36 ORIGINAL PAGE *3 

OF POOR QUALITY 

Thruster Design Impact 

Combining Eqs. 6 and 23 yields the following single equation describing 
the performance of a given thruster design , 


* 



This equation. along with the data just presented suggest, that the performance 

of any thruster design depends on jthe val ues of four physical constants* 

* • 

Gp» C Q , _fg and f c ; and two operating parameters, m and V Q . ...The effect of 
some of these parameters on performance was investigated analytically through 
Eq. 32 by cho.osing the following set of values as the standard set, then 
varying them one at a time to determine their effect on the traditional per- 
formance curves. 


Standard State Parameters 
m - 1000 mAeq 
f B =0.6 

f C = 0 ’ 1 
V D = 50 v 

C Q = 3.0 Aeq" 1 

e* = 50 eV 

Figure 14 shows the effect of- the extracted ion fraction on performance . 
As expected, this parameter has a strong effect on the performance. As seen 


These parameters are constant in the Sense that they are independent... of 


the propellant utilization and flow rate. 
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in this figure, changes .in fg shift the performance curve up or down but do 
not change its shape. Clearly it is desirable to have fg be as large (near 
unity) as possible. 

The effect of the parameter C Q on performance is given in Fig. 15. 

This parameter also has a strong effect on the performance. Indeed, it is 
this parameter which primarily. .determines the shape of the performance curve, 
with larger values of C Q corresponding to curves with more sharply defined 
"knees." From the definition of C Q given in Eq. 24 we see a number of ways 
in which the value of C Q may be. increased. For example C Q may be increased 
by using. a propellant gas characterized by a larger inelastic collision 
cross, section a Q , and a larger atomic mass (resulting in a lower neutral 
velocity v Q ). One may also increase C Q by decreasing the grid transparency 
to neutrals, $ . This must be done, of course, without increasing the- 
accelerator impingement. current. For thruster designs with non-uni foim beam 
profiles tailoring the accelerator grid hole size to match the radial cur- 
rent density profile might be a useful way to minimize <p . Also, three grid 
systems might be expected to have smaller values of <|> than two grid systems. 

Most importantly C Q may be increased by. increasing .the parameter «- e . 

This parameter corresponds to the average distance a primary electron would, 
travel in the discharge chamber before being collected by an anode surface 
provided it had no inelastic collisions. As mentioned earlier the primary, 
function of the magnetic field in all thruster designs is to increase thi s _ 
length. In cusped field thrusters primary electrons are lost to the anode 
through the cusps. Thus sl q may be increased by decreasing the number of 
cusps at anode potential ,. but only up to a point. Reductions in anode cusp 
area below a certain limit will result in unstable operation of the 
discharge.^ 9 
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Equation 24 also suggests that C Q may be increased by masking down the 
area of the grids through which the beam is extracted, A ... However, decrees- 
ing.Ag in this manner will lead to a large reduction in the extracted ion 
fraction, and therefore an overall reduction in performance. 

The effect of propellant flow rate on performance is shown in Fig. 16.. 

In general, higher flow rates produce better performance. The maximum flow 
rate, however, is limited by the ability of the accelerator system to ex-, 
tract the ion current produced (i.e. to extract f Q J p ) . The effect of flow 
rate on performance is less dramatic for thruster designs characterized by 
larger values of C Q as shown in Fig. 17. High values of C Q would therefore 
be particularly, desirable in thrusters designed to be throttled. 

Finally, the effect of c p on performance is illustrated in Fig. 18. 

This parameter merely shifts the performance curves up or. down. The amount 
of the shift increases for smaller valuesof f g . For thrusters that use 
hollow cathodes the efficiency of the cathode operation (characterized by 
V C in Eq. 25) has a strong effect on the value of e*. Further, it is the 
ratio V C /V D that is important. Thus larger values of V Q should produce 
smaller. Ep , however, large discharge voltages are undesirable for a number 
of other reasons. Consequently, a trade off between these conflicting re-„ 
quirements in discharge voltage is necessary. 

Neutral Loss 

Writing Eq, 20 as, 

n Q = m - J B = m (! - n u ) , ( 33 ) 

and substituting this into Eq. 31, and Solving for n Q ' yields, 
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This equation gives the value of the neutral loss rate as a function of the 
beam ion energy cost. For a specified thruster geometry and discharge 
voltage C Q , e p , fg, g c and V Q may be taken to be approximately constant as 
shown earlier. .Thus, since the propellant mass flow rate doesn't appear on 
the right-hand side of Eq. 34, this equation predicts that the neutral loss 
rate n Q is independent of the flow rate for constant values of e This 

same conclusion was reached originally by Kaufman 20 in his constant neutral 
loss rate theory. 


Cone! usions 

A very .simple thruster performance model has been developed which is 
easy to apply as an aid for the design of new thrusters. This model describes 
the performance in terms of four parameters: the plasma ion energy cost 

U p ), the extracted ion fraction (f„), the ion current fraction to cathode 
potential surfaces ( f^.) , and the discharge voltage (VJ. 

The equation developed to describe the behavior of the plasma ion energy 
cost agrees well with a variety of experimental results. The model and ex- 
periments suggest that the direct loss of primary electrons to the anode is 
the primary factor determining the behavior of the plasma ion energy cost 
and the shape of -the performance curves. 

The extracted ion fraction may be ..taken to be Independent of the pro- 

pell ant flow rate and utilization. It does, however, depend on the propellant 
type and discharge voltage. 

For design purposes the thruster performance may be taken to depend 
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only on four physical constants, C Q , e*. f g and f c and two operating 
parameters, m and V g . 

Improved thruster performance should be characterized by high ex- 
tracted Ion fractions, f g , large values of the parameter C Q , small values 
* 

of the parameter e p . 

Future work should focus on determining the dependence of the extracted 
ton faction f g and the primary electron containment length, & e , on the 
thruster geometry and magnetic field configuration. 
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THE EFFECT OF DISCHARGE CHAMBER WALL TEMPERATURE 
ON ION THRUSTER PERFORMANCE 

Introduction 

For many years it has been assumed that neutral atoms in an ion thruster 
discharge chamber move at a velocity determined by the mean discharge chamber 
wall temperature. ’ ’ 6 By applying this assumption it has been possible .to 
make useful predictions of neutral atom density, levels and charge exchange 
ion production rates in ion thrusters. Application of the assumption also 
makes it possible to predict that the performance of a discharge chamber can 
be improved by lowering .the mean discharge chamber wall temperature. The 
extent of the performance improvement that could be induced by such a wall 
temperature reduction has however not been demonstrated. Further, the validity 
of the presumed relationship between discharge chamber wall temperature and 
neutral atom velocity has not been verified experimentally. 

The. intent of the work described herein has been to demonstrate the 
extent of performance improvements induced by reductions in wall temperature 
and to show the discharge chamber model developed in the preceding section 
of this report describes the magnitude of these improvements correctly. 

Various perturbing effects that might cause the atom temperature to differ 
from the mean wall temperature are also examined theoretically. 
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Theory 

The direct measurement of the velocity distribution function or the mean 
velocity of neutral atoms In an Ion thruster discharge chamber is difficult 
to accomplish but Indirect measurements based on discharge chamber performance 
can be used to Infer the neutral atom velocities. In the present case it is 
assumed that the neutral atoms are in thermal equilibrium with discharge 
chamber walls and that the relationship between the neutral atom mean velocity 
and measured discharge performance variables is described in the model pro- 
posed by Brophy. It is argued that agreement between measured changes in 

performance induced by. wall temperature changes and those predicted by the 

model suggest both the. model and the neutral atom/wall thermal equilibriation 
assumption are correct. 

While Brophy's model in its complete form is intended to predict cornier 
discharge chamber performance curves, only a portion of the model will be 
needed here* This portion describes the energy cost of a plasma ion (e ) 
which is the average energy required to produce ions in the plasma. To 

determine this parameter at each operating condition the following relation- 
ship is postulated: 

G P = e p /P o (35) 

where e* is the average energy required to produce an ion in the plasma 
assuming: 1) no primary electrons are lost to the anode and 2) the primary 
electrons have inelastic collisions with ground state neutral propellant 
atoms only. The value of this parameter reflects the facts that energy 
supplied to the plasma by primary electrons is lost as a result Of atomic 
excitation reactions and Maxwellian electron losses to the anode. Typi cally 
G p is a few times the first ionization potential of the propellant atoms so 
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it Is a function of the. pfopel lant. It is also a function of discharge 
voltage (i.e. primary electron energy). 

The parameter P Q in Eq. 35 represents the probability that a primary 
electron will have an inelastic collision with a neutral atom before it is 
lost to the anode. Collisions. between primary electrons and ions are ■ 
neglected because neutral densities are typically an order of magnitude 
greater than the. ion densities. However, at high propeLl ant efficiences 
(high plasma densities) these ionic collisions probably become important and 
the model would be expected to be less accurate there. The probability _P Q is 
determined by recognizing that it is the compliment of the probability that a 
primary electron will escape to the anode before losing its energy. This later 
probability is given by the survival equation, 7 hence 


where o 0 is the total inelastic collision cross for atom/primary electron 
collisions at the .energy of the primary electrons, n Q is the atom density, 
and is the average distance a primary electron would travel in the absence 
of collisions before being lost to the anode. This distance (a ). will be 
referred to herein as the primary electron containment length. 

Considering only, singly charged ions, the neutral flow rate through the 
grids -(n 0 ) is simply the difference between the propellant input flow rate 
(ni) and the ion beam current (Jg) extracted from the thruster 

n o = m - J B (37) 

where it is assumed that each of the quantities in the equation are in units 
of amperes equivalent (A eq.). The_neutral loss rate from the ion source may 
be expressed, however, using the expression for free molecular flow through 
a sharp-edged orifice. 7 
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In this equation n Is the atomic density, e is the electronic charge, A 

w 9 • 

is the area of the extraction grid system exposed to the discharge plasma 
and <j) Q is the effective transparency of the grid system to neutral atoms. 

The neutral atom velocity (v Q ) in Eq. 38 is the mean thermal velocity for 
atoms having a Maxwell -Bol tzmann distribution function 



where k is Boltzmann's constant and m is the mass of a propellant atom. The 
neutral atom temperature. T Q is determined by the state of collisional 
equilibrium achieved between these neutral atoms and the other particles or. 
surfaces with which they interact, in this case, the ions, electrons and 
discharge chamber wall surfaces exposed to the plasma. .Because of the low 
atom and ion densities in an ion source discharge chamber, which are typically 
of order 10 12 cm" and 10 11 cm respectively, atoms will have frequent 
collisions with the walls and one would therefore have the first order 
expectation of an atom temperature near the mean wall temperature (T = T w ). 
Using this equality together with Eqs. 35 to 39 one obtains 


where 
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The suitability of Eqs. 40 and 41 to describe correctly the relationship 
between the energy cost of a plasma ion (e^) and the neutral density 
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parameter m(l - n u ), where n u is the propellant efficiency, has been demonstra- 
ted experimentally by Brophy in tests reported In the preceding section of 
this report. It has been shown that changes in grid transparency ((|> o ), grid 
area (A g ) and propellant type (through the cross section [o 0 ] and the pro- 
pellant atomic mass [m]) are correctly modelled by these equations. 

In this -paper experiments, In which the relationship between energy, cost 
of plasma ions and the neutral density parameter is quantified in terms of 
the coefficient e* and C Q , are described. It is argued that e* is independent 
of wall temperature and that C Q is affected by it only throug.- the explicit 
dependence shown in Eq. 41. Changes in the values of C Q predicted by Eq. 41 
due to changes in wall temperature are compared to measured changes in C 
induced by wall temperature changes in an operating ion source. 

Apparatus and Procedure 

In order to study experimentally the influence of discharge chamber 
wall temperature on discharge chamber performance, the ion source shown 
schematically in Fig. 19 was built. It is similar to the source shown in 
Fig. 2 which was built to study the effects of discharge voltage, propellant, 
grid extraction area and transparency on discharge chamber performance. The 
primary difference between these sources is the cooling coil installed for 
this experiment (Fig, 19). This coil, which was made using 1.3 cm dia. cop- 
per tubing compressed until it deformed plastically onto the steel shell, can 
be connected to a liquid nitrogen source. The liquid nitrogen flow rate 
through the_coil is .adjusted to control the steel shell temperature at values 
ranging from * 80°K to * 500°K. The heat transfer area and contact thermal 
conductance are sufficiently high so wall temperatures measured using the iron- 
constantan thermocouples shown in Fig.. 19 are within * 10°K of the liquid 
nitrogen boiling point (77°K) when the ion source is operated at maximum power. 
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Fig, 19. Liquid Nitrogen Cooled Ion Source Schematic 
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The ion source of Fig. 19 has a 15 cm inside -diameter, steel shell but 
the screen grid is masked down by the magnet support flange so the ion beam 
diameter is only 12 cm. The magnetic field in the discharge chamber is 

produced by the upstream 1.9 cm x 2.6 cm rectangular magnet, the central 

ring magnet assembly and the downstream ring magnet assembly shown in Fig. 19. 

Each of the magnet assemblies shown has been made by placing 1.9 cm x 1.3 cm 

x 0.5 cm thick samarium cobalt magnets end-to-end and side-to-side ..to form 
the rings and the rectangle respectively. The flux density at the surface ... 
of the magnet assemblies is 0.27 T and they are arranged so the inner surface 
of the central ring magnet has a polarity opposite to those at the inner 
surface of the rectangular and downstream magnet assemblies. The central 
ring magnet assembly is insulated from the strip of 0.13 mm thick steel that 
serves as the discharge chamber anode shown in Fig. .19 by a piece of 0.076 mm _ 
thick flexible mica. Both the magnet assemblies and the anode are held in 
place by magnetic attraction forces between the magnets, the steel anode and 
the steel shell. 

In designing and building the discharge chamber special attention was 
given to features that would minimize temperature differences between points 
on the discharge chamber wall and between the walJs and the grid assembly. 

For example, the steel shell, is 0.5 cm thick and joints in the shell are 
silver soldered together. The magnet support flange, which is made of 
0.18 cm thick copper, is also soldered to the shell. A 0.7 mm thick steel 
ring is sandwiched in close mechanical contact between the copper flange and 
the magnet. ring. The screen grid is clamped directly to the copper flange. 

The molybdenum screen grid is 0.43 mm thick, and has a physical transparency 

of 68%. The molybdenum accel grid is 0.53 mm thick and has a physical 
transparency of 57%. The screen and accelerator grids were held at + lkV 
and -0. 5kV respectively during the conduct of the tests. 
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A pair of -0,25 mni dia., 8 cm long tungsten wires, connected in parallel 
and heated to thermionic emission temperatures by passing an alternating 
current through them, serves as the cathode. The maximum cathode power 
applied at the highest electron emission condition required for these tests 
was * 200 w. The source was designed so that either argon or krypton pro- 
pellant could be fed into it through the reverse flow feed line shown in 
Fig. 19. Flow rates of -500, 1000 and 1500 mAeq. were used for both pro- 
pellants. With- argon the discharge voltage was maintained at 50 v and with 
krypton 40 v was used. These. .voltages were selected as sufficient to yield 
stable discharge operation yet low. enough to limit doubly charged ion 
production. 

A special feature of the thruster was that the screen grid, downstream 
ring, and rectangular magnet assemblies, steel shell and magnet support flange 
were all connected in parallel to a power supply that could be used to bias 
these surfaces relative to cathode potential. By biasing these surfaces 
* 30 v native of cathode potential all discharge chamber electrons could 
be repelled from these surfaces and as a result the ion current to them could 
be measured. Adding this current to the beam current one obtains a value for 
the ion production rate (J p ). The error between this measured ion production- 
rate and the true one Is. equal to the ion losses to the anode which have not 
been. included in the measurement. This error has been estimated to be less 
than a few percent based on the effective area for ion loss to the anode 
compared to other ion loss areas. Using the measured ion production rate 
(Jp) the energy, cost of a plasma ion (ep) can be computed using the expression 
derived -in the pr.eced.ing section of this report 

_ ( J D - V V D 


( 42 ) 
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In order to conduct the tests, the ion source was first allowed to 
stabilize thermally at a low discharge current level and a fixed propellant 
flow rate and discharge voltage. The discharge current was then increased 
in increments by increasing the cathode filament current while the flow rate 
and discharge voltage were held constant. At each discharge current the beam 
current and ion current to the steel sheel and screen grid were recorded along 
with the discharge current.. This procedure was repeated at three different 
flow rates for each wall temperature condition for argon and krypton propel- - 
lants. The data were evaluated by computing the total ion production current 
as the sum of the beam current and ion current to the screen grid ancLshell .. 
to find the total ion production rate. Equation 42 was then used to determine 
the energy cost of a plasma ion. A parameter proportional to the neutral 
density was also computed as the difference between the neutral propellant 
flow rate in Aeq and the beam current. 

Error Considerations 

Since .the temperature of propellant atons in the discharge chamber is 
determined by collisions with wall surfaces and influenced by discharge plasma 
ions and electrons there are two general sources of errors that could cause 
the atom temperature to differ from the measured mean_walL temperature. They 
are errors due to 1) variations in wall temperature that complicate the 
computation of the appropriate mean temperature from measured values and .... 

2) perturbations in. the atomic temperature induced either by atomic collisions 
with electrons and ions or an influx of atoms having energies that differ 
markedly from those associated with the mean wall temperature. The magnitude 
of errors that might be expected from these two effects will now be con- 
sidered separately. 
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Wa ll Temperature Variation Effects 

The worst case wall temperature variations are realized when the dis- 
charge chamber and cathode are being operated at maximum power. At this 
operating condition it is estimated that 200 watts of power could be radiated 
from the cathode and the plasma. This power would generally be assumed to 
be distributed uniformly over the interior surfaces of the chamber. Other 
power losses to the walls are associatecLwith electron and ion losses to these 
surfaces. Electrons are assumed conservatively to have an average energy of 
10 eV and to flow to the anode at a rate given by the discharge current 
(4 A maximum). Ions on the other hand are assumed to go primarily into the 
beam, the grid webbing and the grid support flange. 

A consideration of the thermal conduction characteristics of the various.-., 
component parts of the discharge chamber along with the heat loads on these 
parts suggests that the screen grid, anode and cathode surfaces can be ex- 
pected to perturb the mean wall temperature away from the values measured on 
the steel shell.. The worst case temperature profiles over each of the surfaces 
have therefore been determined and the impact of these temperatures on the 
mean wall surface temperature. has been evaluated. The screen grid has the 
greatest affect because 1) it has a low thermal conductivity as a result of 
its thin, perforated configuration, 2) it has a substantial area and 3) it _ 
has a high distributed heat load as a result of direct ion impact over its 
surface. 

In modelling the grids thermally it was assumed that the extracted beam 
current was at its maximum value of 0.9 A. Since the extracted ion fraction 
for this ion source was typically 60% of the total ion production this means 
0.6 A of ion current goes to interior surfaces of the discharge chamber. 
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Assuming a screen grid. transparency to. ions of 80% gives an ion current of 
0.22 A- to the screen grid webbing, and the remainder 0.38-A to other surfaces 
accepting ions. It was also, assumed that the ion current density to the 
screen grid webbing was the same as that through the screen grid .holes and 
that ions striking the webbing had been,accelerated from a worst case plasma 

potential 50 v) to screen grid bias potential (- 30 v). This heat load 

combined with the radiated one induced a worst case temperature at the screen 
grid centerline that was ^ 125°K above the temperature at the inner edge of 
the magnet support flange (Fig. 19). It was next assumed that all ions pro- 
duced in the discharge chamber that did not go into the beam or. onto the screen 
grid webbing (0.38 A) were directed uniformly onto the magnet support flange 

at an energy of 80 eV. Because this plate was made of copper and was thicker 

than the screen grid, the worst case temperature difference computed between 
its inside edge (6 cm radius) and its outer one (7.5 cm radius) was only ^ 7°K. 

The difference in temperature between the anode surface and. the steel 
shell was computed by considering the heat flux to the interior surface of 
the anode resulting from both the electron current and radiated power from 
the cathode and plasma. .It was assumed that the relatively soft mica was 
pressed sufficiently firmly between the magnets and the anode itself so 
contact thermal resistances were small and the major resistance to heat 
transfer occurred in the mica. Under these conditions the temperature dif- 
ference between the anode and the steel shell was approximately 1°K. The 
cathode wire temperatures were assumed to be the value required to assure 
thermionic emission (-v 1475°K). 

When the area weighted mean. temperature -of the surfaces exposed to the 
plasma was computed from the worst. case component temperatures just cited 
it was found to be only ...> 7®K above the mean steel shell temperature. Hence 
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under worst case thermal loading conditions the mean wall temperature would 
be expected to be 7°K above the mean steel shell temperature. During the 
actual conduct of the test the two. thermocouples shown in Rig. 19 agreed .to 
within about 5°K and they were found to exhibit worst case drifts of ^ 10°K 
during the conduct of a test in the temperature range 80 to 160°K and ^ 5*K 
during the conduct of the higher temperature tests. -Based on these measure- 
ments and calculations it is suggested that mean wall temperatures agree with 
measured ones to within t 10°K for low temperature tests and t 5°K for higher 
temperature ones. 

Colli signal Effects 

The neutral atom temperature could conceivably differ from the mean 
wall, temperature because of collisions with electrons and ions which have 
energies that differ substantially from those associated with the unperturbed 
atom temperature. Electrons, because of their mass mismatch with the atoms, 
would not however be expected to perturb the atomic temperature significantly 
even though their temperature is typically quite high 10 eV = 110,000°K). 

Ions on the other hand could affect the atom temperature because they have a 
larger elastic, energy transfer cross-section, because they can undergo charge . 
exchange reactions with atoms and because they can be accelerated to substan- 
tial energies by electric fields in the plasma. For typical discharge chamber 

configurations the ion energies can generally reach the electron temperature 
~ 22 

(< 10 eV) as determined by the Bohm criterion but will probably not exceed 
it. Typical elastic collision cross sections between argon or krypton ions 

6 o - 23 

and atoms lie below 10 A at the low (<...10 eV) energies involved. For a 
typical ion density of 10 11 cm" 3 this corresponds to a mean free path greater 
than 100 m which is 700 times the discharge chamber dimensions. Hence collisions 
with vial 1 s should dominate ’nd elastic collision effects should be negligible. 
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For the .charge exchange process 2 **’ 25 the argon and krypton cross sections 

°2 

are higher (of the order of 50 A ) but this still implies mean free paths 
that are of the order of 100 times the discharge chamber dimensions so this 
colli slonal event would still be unlikely compared to wall collisions. 

Finally.it is noted that ions that go to discharge chamber walls rather 
than being extracted into the beam are accelerated into the walls through a 
plasma -to -wall potential difference which is of the order of several tens 
of volts. These ions while they could recombine and come into thermal 
equilibrium with the wall before they leave it might also be expected to re- 
combine and recoil from it with some fraction of their incident energy. While 
it was hoped that the error due to recombined ions recoiling from the walls 
with energies above those expected for particles in thermal equilibrium with 
the walls would be small, no conclusive evidence was found in the literature 
to substantiate a supposition that this would or would not be the case. 

Results 

When the ion source was operated on argon at a wall temperature of 390°K 

and a discharge voltage of 50 v the plasma ion energy cost and neutral density 

parameter data pairs indicated by the data symbols of Fig. 20 were obtained.. 

Through trial and error selection of the parameters c* and C Q ,Eq. 40 was fit 

to these data. The fit, shown by the solid line of Fig. 20, was obtained for. 

C = 3.0 A eq- 1 and e *56eV. The parameter e is determined to first order 
o p P 

by the propellant and the energy of the primary electrons (i.e. the discharge 
voltage) so one would not expect it to be affected by changes in wall tempera- 
ture. Changes in wall temperature should however cause the parameter C Q to 
change in accordance with Eq. 41, and this equation could be applied directly 
if values for all of the quantities appearing in it were known. While most 
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of the values, are known, we unfortunately do not have, a model for the primary 
electron containment length One can still see if the model describes the 
physical situation adequately however by determining the value of C Q appro- 
priate to one wall temperature, then changing the wall temperature while 
holding all else constant and determining if the prediced change in C Q is 
consistent with Eq. 41. This has been done here by taking the valueof C Q 
obtained by curve fitting the data of Fig. 20 (C Q re ^) and computing new 
values of C Q appropriate to any new wall temperature from the following 
equation which has been obtained from Eq. 41. 


C 


o 


C o ref 



(43) 


Doing this one predicts the parameter C Q should be 3.7 A eq." 1 for the case 
where the mean wall temperature is 255°K. Using this value the solid line 
shown on Fig. 21 is predicted. When the ion source was operated at this wall 
temperature the data points shown on Fig. 21 were measured. The agreement 
between these measured data points and the predicted curve is observed to be 


excel! ant. Applying Eq. 43 for a mean wall temperature of 90°K one obtains 
C Q = 6.2 A eq" 1 . This corresponds to the solid curve shown in Fig. 22. When 
plasma ion energy cost/neutral .density parameter data pairs were measured with 
the ion source operating at 90°K the results shown by the. data points on Fig. 22 
were obtained. Again the agreement between theory (solid line) and experiment 
(data points) appears to be good. 

The ion source was also operated on krypton propellant at a mean wall 
temperature of 390°K. In order, to predict the performance with this new 
propellant the constant C 0 was computed from the expression 


C = C r 
o o ref 



( 44 ) 
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where the subscripts Kr and Appertain to argon and krypton respectively. 
Using the appropriate atomic masses and total Inelastic collision cross 
sections for argon at 50 vand krypton at 40 v as obtained from deHeer 14 In 
Eq. 44 one finds C Q * 5,5 A eq' 1 should be appropriate. The baseline 
energy cost per plasma ion e* would -also be expected to change with the 
propellant and discharge voltage. Its value was selected in the present 
case for krypton from experimental data as the plasma ion energy cost at a 
high neutral density parameter i.e. high flow rate/low beam current operating 
condition and was found to be 48 eV. Using C Q = 5.5 A eq' 1 and c * = 48 eV 
in Eq. 40, the solid line shown in Fig. 23 was predicted for operation on 
krypton at 390°K. The discharge chamber performance data measured with 
krypton at 390°K are shown by the data symbols on Fig. 23. Once again the 
agreement between the predicted and measured performance is good. For a 
krypton fed discharge chamber cooled to 255°K and,155°K values of the 
parameter C Q equal to 6.8 and 8.7 A eq’* respectively are predicted by Eqs. 42 
and 43. Using these values as input to Eq. 40 the solid lines of Figs. 24 
and 25 are predicted. Measured data for these temperatures are again shown 
by the data symbols in these two figures. As in all previous cases the 
measured data agree well with the predicted curves. It is noted that dis- 
charge chamber wall temperatures could not be reduced to * 90°K when krypton 
propellant was being used because. the propellant tended to condense. 

The results shown in Figs. 21 through 25. all show sufficiently close 
agreement between predicted and measured performance results to suggest the 
effect of wall temperature is properly modelled In Eqs. 40 and 41. Changes 
in the constant C Q of +10% were generally sufficient to cause predicted 
curves in the figures to depart from tr.ameasured data to an easily noticeable 
degree. As a result it is suggested that the parameter C Q is accurate to within 
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about t 10% and the wall temperature effect is modeled accurately by Eq. 41 
to within t 20%. 

The fraction of the ions produced in the discharge chamber. that were 
extracted into the ion beam was also measured at each operating condition 
for these tests. These extracted ion fractions were relatively constant 
over the full range of neutral density parameters investigated taking on 
values essentially the same as those shown in Figs. 12 and 13. The ex- 
tracted ion fractions were independent of discharge chamber wall temperature.. 

The results presented thus far, while they suggest the validity of the 
model proposed, do not show the extent to which the more traditional perfor- 
mance curves of energy cost per beam ion vs. propellant utilization are al- 
tered by changes in discharge chamber wall temperature. A typical example 
of the magnitude of the influence of .wall cooling is shown in Fig. 26 and 
it is seen to be substantial. A point which should be made about the results 
of Fig. 26 is that the beam ion energy costs could be lowered substantially 
below those observed in the figure by using small hole accelerator grid 
optics and by increasing the propellant flow rate above the value used to 
obtain the data shown. Finally, it is also noted that the change in the 
traditional performance curves predicted by the complete theoretical model 
described in the previous section of this report is essentially identical to 
that shown by the experimental data in Fig. 26. This occurs because the 
high degree of agreement realized in the data of Fig. 20 through 25. 


Conclusions 

Substantial improvements in discharge chamber performance are induced 
by reductions in the mean discharge wall temperature. The extent of these 
improvements are predicted by the simple discharge chamber model described 
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by Eqs. 40 and 41. Based on the high degree of agreement between the pre 
dictions of this model and experimental results it is concluded that the 
assumptions of the model are valid. These assumptions are 1) that the 
atoms present in a typical ion thruster discharge chamber move at a mean 
thermal velocity determined to first order by the mean discharge chamber 
wall temperature and 2) that this mean velocity is not perturbed signifi 
cantly by collisional or recombined-ion recoil effects. 


LOW MAGNETIC FIELD STRENGTH DISCHARGE CHAMBERS 
J. R. Brophy 

The discharge chamber model discussed in the preceding sections of 
this report was developed for cusp-type discharge chambers where electrons 
are lost predominantly into the cusps along magnetic lines-of-force. For 
discharge chamber configurations in which electrons must diffuse across 
magnetic field lines the field strengths can be considerably less than those 
found in the cusp-type chambers. Indeed, they must be considerably weaker 
in the diffusion-type discharge chamber in order to effect proper operation. 

In spite of these differences many of the operational features of the two 
discharge chamber types are similar. In both cases primary electrons are 
injected into the discharge chamber plasma, ions are created by electron 
bombardment, a fraction of the ions produced are extracted into the beam, 
and Maxwellian electrons are collected by the anode. In the preceding 
sections of this report it has been demonstrated that the performance of a 
cusp-type thruster is strongly dependent on the loss of primary_electrons 
through the cusps to the anode. The similarity in the shape of the perfor- 
mance curves for both the cusp-type and diffusion- type discharge chambers 
suggests that the loss of high energy or primary electrons. is also occurring 
in the low field strength diffusion-type configurations. It can also be 
argued that this is reasonable from an energy conservation standpoint. 

Consider the situation at high propellant utilizations where increasing the 
cathode emission current produces little additional beam current. This 
suggests that few additional ions are created. Yet, the energy contained 
in these additional primary electrons must go somewhere. This energy may 
go into excitation collisions with ions, increasing the Maxwellian electron 
temperature or it may be carried directly to the anode by the primary electrons 
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It is this third possibility that seems the most likely based on the results 
of this report. In fact, calculations by Longhurst 26 support this conclusion. 
For this hypothesis to be viable, however, a. model describing primary elec- 
tron transport across magnetic field lines without energy degradation must 
be identified and verified experimentally,. 


EFFECT OF SCREEN GRID POTENTIAL ON PERVEANCE 
J. R. Brophy 

Introduction 

27 

It has been observed that the impingement-limited perveance of a two- 
grid accelerator system is a function of the ratio of discharge-to-total 
voltage (Vg/Vj),_ provided the screen grid is at cathode potential and the 
plasma is. near anode potential. These observations indicate that the per- 
veance decreases as the ratio Vg/V-j- is increased. It is believed that this 
effect is predominately the result of the potential difference between the 
plasma and the screen grid (aV), which in the above case is .approximately 
equal to Vp. Other effects, such as. a change in Maxwellian electron tem- 
perature or primary electron energy, resulting from a change in discharge 
voltage .are not believed to bo important. 

To investigate the effect of aV/Vj on perveance experimentally a 15 cm. 
dia. ring cusp ion thruster was configured with the screen grid electrically 
isolated from the thruster body. This allowed the screen grid to be biased 
negative of cathode potential while the discharge voltage (V D ) and total 
voltage (V T ) were held constant. In this case the potential difference 
between the screen grid and the plasma is given approximately by, 

AV = V D + | V s | , (45) 

where both voltages are measured relative to cathode potential and |V^| is 
the magnitude of the negative bias applied to the screen .grid. Earlier 
attempts to measure the effect of the screen grid-plasma potential difference 
on perveance by varying the discharge voltage while holding the total volt- 
age, constant have been limited by the range over which the discharge voltage 
Could be varied. By biasing the screen grid in the manner described above. 
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the impingement-limited perv.e.ance can be measured over a large range of 
AV/16,- while the discharge voltage, total voltage and net-to-total voltage 
ratio are held constant* 

Apparatus and Procedure 

These experiments were conducted on a 15 cm dia, ring cusp ion thruster 
equipped with a set of dished, small hole accelerator grid .(SHAG) optics. 

The physical open area fractions of the screen grid and accelerator grid 
were 0.68 and 0.30, respectively, and the cold grid spacing was .75 mm. 

The following procedure was implemented for the collection of data in 
this investigation. The flow rate of the argon prop! ell ant used was set 
and maintained at 1000 mAeq,and the discharge vol tage was held constant at 

50 v. Total voltages of 500, 1000 and 1500v were investigated at a net-to- 

total voltage ratio of 0.6, At each value of the negative Dias applied to 
the screen grid, the discharge current, was adjusted until the impingement- 
limited perveance condition was realized. This condition was said to exist 
when the impingement current started to rise dramatically. The use of non- 
regulated high voltage power supplies necessitated the use of this somewhat 
imprecise definition of the impingement limit. With these power supplies 
an increase in impingement current causes .increased power supply loading 
which then causes the accelerator grid potential to become less negative. 

This in turn induced a further increase in impingement current and these 
effects cascade. The perveance level at which this process begins is 
designated at the impingement limited perveance. 

Results and Discussion 

The results of these tests are shown in Fig. 27. The parameter J n /\4 3 / 2 , 

B T 9 

where J Q is the beam current at the impingement limit, is proportional to 
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Fig. 27. Effect of Plasma-to-Screen Grid Potential Difference 
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the Impingement-limited perveance. Clearly this parameter Is a decreasing, 
function of aV/V t , and is approximately Independent of the total voltage, .. 

These results, obtained with a 15 cm dla, ion source, .are consistant with 
those found by Rovang 1 on a much smaller test facility. The experiments of 
Ref. 1 were performed on an 8 cm dla. source masked to produce an Ion beam 
approximately 1 cm In diameter centered on the thruster axis. This insured 
ion extraction from a plasma of near uniform density across the active grid 
area. In the present case, however, the plasma density is almost certainly 
non-uniform across the grid area. Thus it is considered significant that 
the results of the present test agree with those obtained by RovangJ 

Since the results of Fig, 27 indicate that the impingement-limited 
perveance decreases as the ratio AV/Vy is increased, it is desirable that 
AV/Vy be made as small as possible. It should be noted, however, that if 
aV is made smaller by biasing the screen grid more positive than the floating 
potential of the plasma a discharge chamber performance penality in the form 
of higher plasma and beam ion energy costs will be incurred. These higher 
ion energy costs will result from the increased removal of energetic electrons 
from the plasma to the screen grid. 

For missions characterized by low values of specific impulse (low net 
accelerating voltages), the requirement. for low AV/Vy might be expected to 
impact the selection of the desired net-to-total accelerating voltage ratio (R), 
Greater flexibility in the selection of this ratio would be afforded if 
AV/Vy could be reduced by decreasing a V. For direct current discharge cham-. 
bers a reasonable limit on AV would be expected to be the discharge voltage/ 
floating potential difference for the reasons mentioned above. For RF dis- 
charge chambers, however, a V could presumably be made as small as a few volts. 


i, 







t 

\ 


I 


I 


I 


.iMLLfcJI 





76 


XENON AND ARGON HOLLOW CATHODE RESEARCH 
Dan Siegfried 

Introduction 

Recent -research has provided a better understanding of the basic 
physical processes underlaying operation of orificed hollow cathodes* 8 and 
has led to the development of an analytical model describing these processes. 
That work was done for cathodes operating on mercury, the propellant. cur- 
rently used in flight-certified ion thrusters. The development of advanced 
ion thrusters which operate on. xenon and argon has resulted in the. need for 
main and neutralizer. hollow cathodes which operate efficiently and and re- . 
liably on these gases. At present, hollow cathodes operating on these gases 
under neutralizer conditions have relatively high electrical power and pro- 
pellant flow, rate requirements. In order to develop a more efficient hollow 
cathode for this application, it would be useful to understand how the use of 
argon and xenon affect the operation of the cathode. 

The objective of the research described here is to determine what effect 
operation with these propellants has on the basic physical processes previously 
identified. for mercury, orificed hollow cathodes and to determine if the 
analytical model developed for mercury cathodes can also be used with argon 
and xenon. In order to do this, experiments similar to those conducted with 
mercury were performed using a special quartz-body test cathode and mounting 
fixture. This apparatus allows the Simultaneous measurement of cathode 
internal pressure, the insert temperature profile, internal plasma properties, 

and the emission currents from various cathode Surfaces, The results of 

theSe experiments will be presented here and compared with calculated. results 
based on the analytical model. 
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Apparatus and Procedure 

The cathode test configuration used In the experiments Is shown In 
Fig. 28. The cathode outer tube was made of quartz to allow direct visual 
observation of the Insert for temperature measurement and to provide 
electrical isolation of the insert for current measurement. The cathode 
used a single layer tantalum foil insert dip-coated in. the chemical R-500.* 
The Insert was 15 to 20. mm long and had a diameter of 3.8 mm; the cathode 
was operated with orifice plates having 0.51 and 0.76 mm diameter orifices. 

A movable Langmuir probe mounted on the cathode axis was used to measure __ 
plasma properties within the cathode. A manometer having diffusion pump 
oil as a workina fluid was used for sensing the upstream stagnation pressure 

in the cathode and the insert temperatures were measured with a micro-optical 
pyrometer. 

The apparatus and procedures used in the experiment were essentially 
the same as those used in the earlier experiments with mercury. The reader 
is referred to the report on those earlier mercury experiments 28 for details 

of construction and procedure, since only the few notable differences will be 
pointed out here. 

Because the recent experiments were run with a non-condensing gas, the 
heater on the manometer was not. needed and war, not used. In addition, gas 
flow rate adjustment was accomplished with a micro-adjustable needle valve 
in the propellant feed line rather than the .vaporizer heater required for 
mercury. Adjustment of the flow rate was, therefore, much simpler than with 
mercury and the throttle valve shown in. Fig. 28 was not necessary 6 This valve 


A double carbonate mixture - (Ba/Sr) C0 3 manufactured by the J.R. Baker 
Chemical Co., Phillipsburg, New Jersey, 
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was left closed In the present tests. 

In the earlier tests, the electrode used in the probe -was a tungsten 
sphere with a diameter of * 0.76 mm. formed, on the end of a 0.25 mm diameter 
wire. Using that probe, perturbation Of the Internal discharge appeared 
to be somewhat greater with argon and xenon than that experienced with 
mercury. For this reason, probes having a smaller electrode surface area 
were used in the present experiment. Spherical probes ranging in diameter 

from 0.48 to 0.58 mm and a planar probe with a surface 0.25 mm in diameter. 
were used. 

Depletion/contamination of the insert was a significant problem in 
the argon and xenon experiments presented here. Obtaining and maintaining an 
insert surface having a relatively low, uniform and stable work function proved 

to be extremely difficult. . This was somewhat of a problem, though much less 

severe,, when the experiments were performed with mercury. In the experiments 
with all of the gases including mercury, the materials and construction, were 
the same; and a great deal of care was taken in fabricating and conditioning 
the insert. However, because of the severity of the problem in the present 
case, two additional precautions were taken which were not used in the mer- 
cury experiments. One was the installation of a liquid nitrogen cold finger 
in .the bell jar near the cathode. The purpose of this was to help remove 
residual water vapor and other condensable contaminants present in the bell 
jar. The cold finger was operated during insert conditioning and during 
cathode operation. In addition, greater care was taken in constructing and 
outgassing the Langmuir probe used in the tests. Although both of these 
procedures seemed to improve the problem of insert degradation to a small 
degree, it now appears that the degradation may be directly related to the 
use of xenon and argon as propellants. This will be discussed in more detail 
with the presentation of the results. 
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The problems with insert degradation mentioned above .affected the 
method in which data was taken during the experiments. The insert work 
function has a significant effect on such cathode operating parameters as 
insert temperature and plasma, properties. Since insert work function tended 
to change during the course of an experiment, an attempt was made to take 
data as quickly as possible to minimize the effect of this work function 
variation., In addition, complete data, sets (plasma property profiles, insert 
temperature profiles, and current distributions) for each operating condition 
were not always recorded simultaneously. The reason for this was related to 
the fact that the effect of a parameter such as flow rate on, say, insert 
temperature could be masked easily by changes in insert condition during the 
time required to record a complete data set. Because of this, there were a 
number of instances in which a. single measurement such, as plasma density or 
insert temperature were made as a particular operating parameter such as flow 
rapidly varied. Such~a procedure has an obvious advantage in showing 
the effect of one parameter on another, but does not solve the problem of 
obtaining complete and consistent sets of- data over a range of operating 
conditions .in the presence of changes in work function. There was little 
that could be done about the latter problem. The considerable scatter found 
in some of the results is believed to be due in a large part to this 
difficulty. 


The Phenomenological Model 

The earlier work done with mercury resulted in the development of 
phenomenological model describing the important physical processes under- 
lying the operation (ii the cathode. One of . the purposes of the present 
work is to determine whether that model applies to cathodes operating on 
argon and xenon. R wi’l be useful, therefore, to begin by stating the 
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key aspects of the model developed for mercury. This will then be used 
as a background for dicussing the results obtained with xenon and argon. 

The experiments with mercury showed that the electrons which exit 
through. the cathode orifice are produced within the cathode both by sur- 
face emission 70%) and by volume ionization processes ^ 30%), The 
surface emission is due primarily to field-enhanced thermionic emission 
from a relatively well defined band on the downstream end of the insert; 
while volume ionization occurs primarily within the region circumscribed 
by the.emitting portion of the insert. Acceleration of surface emitted 
electrons across the plasma sheath results in the formation of a nearly, 
monoenergetic population of primary electrons. The_primaries are essential 
for sustaining the volume ionization process. They represent the sole energy... 
input into the "ion production region": and are estimated to be directly 
responsible for v 65% of the ion/electron pairs produced in the plasma. 

The rest of the ionization is attributable _to the high energy tail of the 
Maxwellian electron population. Since for mercury the cathode internal 
plasma potential is typically ^ 9 v, the ionization is predominantly a 
multi-step process, relying heavily on the production of ions, from inter- 
mediate metastable and resonance states. The ions produced in the "ion 
production region" diffuse out of it at the Bohm velocity. and are accelerated 
across the plasma sheath striking the insert with sufficient energy to heat 
it to the emission temperature. These ions are neutralized at.the insert 
surface and thus complete the current path between the cathode surface 
and the volume-produced electrons which exit through the cathode orifice. 

The processes described above are discussed in detail in Ref. 28 where 
it is shown that they can be represented analytically in a rather simple 
form if two assumptions are made. First, the emitting length of the. insert 
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and, therefore, the length. of the Ion production region is assumed to be on _ 
the order of the energy exchange mean free path for a primary electron. The 
rationale here is that the primary electrons are essential for the processes 
in the region and that as soon as they .sustain a significant energy loss in 
a collision they are quickly thermalized and cease to exist as primary 
electrons. Having defined the length of the ion production region based on 
the primary electron energy exchange mean free path, it is next assumed that 
the plasma properties (plasma density, plasma potential, and electron 
temperature) and the insert surface temperature are uniform throughout this 
region. In addition to. the two assumptions above, an experimentally de- 
termined value is assumed for the Maxwellian electron temperature. The 
reason for this is that there is a priori no simple. means of predicting the 
electron temperature and the model is relatively insensitive to this param- 
eter over its normal range of values. 

Results 

Results will be. presented in four sections in the following order. 

First, typical results will be presented which show the qualitative effect 
of the propellant type and the insert surface work function on cathode 
parameters such insert temperature, and plasma properties. These results 
will provide general information regarding the effect. of propellant type on 
the basic physical processes occurring within the hollow cathode. In the 
next section, inelastic cross-section data will be presented for argon and 
xenon and will be used to shed light on the nature of the collisional 
processes taking place in the ion production region. In the third section, 
detailed results for argon and xenon will be presented which show the 
effect of discharge current and flow rate on plasma properties and insert 
temperature; and a comparison will be made with predictions of the analytical 


model. The final section will present a correlation of press ure-f low rate 
data for argon and xenon. The results of that section can be used to pre- 
dict the cathode internal pressure when applying the analytical model to 
situations where the pressure is unknown, 

Propella nt and Insert Surface Work Function Effects 

The effects of propellant gas and insert surface work function On 
the Internal plasma density profile are shown in Fig. 29 for Similar cathodes 
operating at discharge currents of 2.3 A and an internal pressure .of ^ 4 Torr. 
The solid curves represent results obtained in quick succession with the 
same cathode by simply turning off the xenon supply and turning on the argon 
supply (the two gas supplies were connected in parallel). Since the probe 
and insert surface conditions are believed to be nearly the same in both 
cases,, this provides a direct comparison of the effect on plasma density of 
varying the propellant from xenon to argon. The results show that the peak 
plasma density for xenon is nearly double that for argon and suggests, that 
the ion production region is on the order of 5 mm long for xenon while for 
argon it. is closer to 7 mm. A comparison can also be made between operation 
with mercury (dashed curve) and with xenon (center-dash curve) for similar 
cathode conditions where both cathodes are operating with inserts having 
surface work functions of ^,1.9 eV„ The xenon cathode achieves a peak 
plasma density ^ 505S higher than the mercury and the ion production region 

for the xenon appears to be about twice the 1 to 2 mm length found with 
mercury. 

The effect of insert surface work function on plasma density can be 
seen by Comparing the difference between the two curves for xenon in Fig. 29. 
The main consequence of a lower work function appears to be a reduction in 
the length of the emission/ion production region. The difference in peak 
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plasma density between~the two xenon curves in Fig. 29 is not believed to 
be significant. It is less than typical scatter found in the data and, 
in fact, other results have indicated that the peak plasma density may 
increase with increases in surface work function. 

A similar comparison is presented in Fig, 30 for insert temperature 
profiles, again with propellant type and surface work function as parameters. 
The results of Fig. 30. suggest that the peak insert temperature is relatively 
independent of the propellant type for a given value of the insert surface 
work function. It is also worth pointing out in Fig. 30 that the insert 
temperature peaks correspond approximately with the plasma density peaks 
in Fig. 29 and that the emission lengths indicated by the temperature profiles 
are also consistent with the length of the ion production region indicated 
by the plasma density profiles in Fig. 29. 

Figure 31 shows the effect of propellant type and insert work function 
on plasma potentials inside the cathode. The difference between the two 
xenon curves indicates that the plasma potential increased significantly with 
surface work function. The two solid curves shown in Fig. 31 are again the 
results obtained using the same cathode and simply switching from xenon to 
argon. A comparison of. these two curves shows that under similar surface 
work function conditions (<f> $ 'v 2.5 eV) that the internal plasma potential 
with argon is ^ 3 volts greater than that with xenon. This is consistent 
with the 3 volt difference both in first excitation and first ionization 
energies between the two gases. On the other hand, the xenon and mercury 
cathodes each operating under similar conditions with a similar surface 
work function ($ s 'v 1.9 eV), the-.internal plasma potentials are similar 
(Fig. 31, dashed and center-dashed curves) while on the basis of either 
first excitation or ionization potential xenon might be expected to exhibit 
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a higher plasma potential than mercury. It is unclear as to why this is 
not the case. 

The effect of insert surface work function and propellant type on the 
electron temperature is shown in Fig. 32. In the upper plot of Fig. 32 the 
solid curves show the electron temperature of xenon plotted as a function 
of discharge. current with surface work function .as a parameter, while the 
lower plot shows similar results for xenon plotted as a function of internal 
pressure. The dashed curves in Fig. 32 are for a similar cathode operating 
on mercury and having an insert with a work function of ^ 1.9 eV. No curve 
is shown for the argon data because of the few data points available. There 
are two significant points, to be noted in Fig. 32: 1) the electron temperature 

for xenon and mercury are similar with average values near -v 0.8 eV, when the 
surface work function is vl.9 eV, while argon has a significant! y higher, 
electron temperature near 1.3 eV; and 2) the electron temperature increases. 
as the insert surface work function increases (indicated by the. xenon curves). .. 

In addition, to the propellant effects noted in Figs. 29 to 32 there are 
some qualitative differences which were observed during the experiments. and 
are worth noting here. One of these is the apparently greater sensitivity, 
of the insert to depletion or contamination seen when operating the cathode 
on xenon and argon as compared to operation with mercury. The reason for 
this is not clear. One possibility is increased surface camage/depletion 
due to sputtering by ions returning to the surface from a plasma at a higher 
plasma potential, an explanation that appears particularly likely in the case 
of argon. Another possibility is increased sputtering due to higher ion 
current densities to the surface, particularly with xenon. The ion flux to 
the surface is directly proportional to the local plasma density and 
inversely proportional to the square root Of the atomic mass. Xenon 
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exhibits higher plasma densities than does mercury under similar conditions 
and has a lower atomic mass, both of which result In greater Ion fluxes to 
the surface than those seen with mercury. These higher ion fluxes to the. 
surface could also be a factor in enhancing the depletion rate for the low 
work function oxide on the insert. 

Another difference seen when operating with argon and xenon was a 
greater tendency for the electron emission region to locate itself at some 
upstream point on the insert. Although, as with mercury, the most stable 
emission site is apparently at or near the downstream end of .the insert, a 
relatively large number of. instances of emission from other sites was ob- 
served in the case of argon and xenon. In addition, when electron emission 
was established at an upstream location with argon and xenon, the alternate 
site was relatively more stable and less likely to shift back to its normal 
downstream location than it was. in similar situations observed with mercury. 
Emission from regions other than the downstream end of the insert may be a 
reflection of the greater insert depletion rate discussed above, because 
such behavior was usually observed during operation under higher surface 
work function conditions. 

Another related phenomenon observed with argon and xenon was the tendency 
to form localized emission spots one to two millimeters in diameter. Although 
these were also observed occasionally during operation with mercury, their 
occurrence was more common with argon and xenon. This was in part due to 
the fact that a considerable amount of the operating time particularly in 
the initial xenon experiments was at a relatively high internal pressure? and 
the occurrence. of these localized spots exhibited a definite correlation 
with operation at high internal pressure. In fact, initiating such 
localized emission could b e readily accomplished by operation at pressures 


I 


91 

above * 10 To nr. The results -of the experiments suggest that operation at 
such, high pressures should probably be avoided because. the Intense local 
activity accelerates depletion of low work function materials. from the Insert. 

The results presented In Figs. 29 to 32 and discussed above suggest that 
the basic physical processes for cathodes operating on argon and xenon are 
probably very similar to the ones found in mercury cathodes. That is surface 
emission normally takes place from a relatively well defined region near the 
downstream end of the insert, predominantly as a. result of field-enhanced 
thermionic emission. The emission from this region coincides with and, 
indeed, sustains the intense excitation and ionization in the volume 
adjacent to the surface emission region. The length of the emission/ion 
production region is different for the three propellant types as are the 
magnitudes of the plasma properties but, at least at first glance, the basic, 
physical processes would appear to be similar. 

Colli siona l Processes in the Ion Production Region 

In the next section more detailed results will be presented which 
show the effect of discharge current and mass flow rate on the operation 
of the hollow cathode with xenon and argon. Those results will provide 
an opportunity for comparison with the predictions of the analytical model. 
However, in order to make calculations with the model, it is first necessary 
to know the total inelastic cross-sections describing electron/atom and 
electron/ion collisions at energies near the primary electron energy. This 
energy is numerically the same as the plasma potential. These cross-sections 
are needed to determine the primary electron energy exchange mean free path 
and also to provide general information about important collisionai processes 
taking place within the cathode. This section will present the necessary 

cross-sections for argon and xenon and discuss them in terms of the relevant 
collisionai processes. 

c -a- 
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The total Inelastic cross-section for-electron impact is shown as a__ 
function of. electron energy in Fig. 33 for argcn and xenon. For energies 
below the first ionization level (15.7 eV for argon; 12.1 eV for xenon,) 
the data come, exclusively from Ref. 29. They are .presented in the form of 
a composite excitation cross-section containing all of the excitation levels 
lumped together. For energies above the first ionization energy, the 
ionization cross-sections (Ref. 30) are added to the excitation cross-sections 
(Ref. 29) to obtain the results shown in Fig. 33. 

Argon and xenon cathodes both behave very much like cathodes, operating 
on mercury. ...However, in the case of xenon there appear to be significant 
differences in the details of the important colli sional processes. For 
example, Fig. 34 shows the inelastic mean free path for a primary electron 
colliding with xenon atoms as a function of xenon atom density and with 
primary electron energy as a parameter. The same plot for mercury is shown 
for comparison in Fig. 35. For typical operating conditions (neutral 
density = 3 x 10 16 cm" 3 and plasma potential 'v 9 v) xenon is seen to have 
an inelastic mean free path about 25 times the 'v 1 mm value. shown for mercury. 
In addition. Fig. .34 shows that for xenon the mean free path exhibits a very 
strong energy dependence particularly for electron energies in the range 
observed ir a typical cathode (8.5 - 10 eV). For xenon, the mean free path 
is very sensitive to energy over this range because the energies are near 
the excitation threshold energy, while for mercury electron energies in the 
8.5 to 10 eV range are considerably above the first excitation level (4.8 eV). 

In the model developed for jnercury cathodes,, the energy exchange mean 
free path for a primary electron is used as a criterion f or. J:he length of 
the emission region. This quantity depends not.only on the inelastic mean 
free path (Xj^ shown in Figs. 34 and 35), but also on the energy exchange 


ELECTRON ENERGY (eV) 


Fig. 33. Argon and Xenon Inelastic Cross. Secti 





PRIMARY ELECTRON ENERGY (eV) 



Electron Inelastic Mean Free Paths in Xenon 





96 


mean free path associated with elastic scattering CQlllsions..between primary 
electrons and the low energy Maxwellian electrons U e ); specifically 

V = (1/x In + 1/A e )-1 * (47) 

The electron scattering, mean, free path (a ) is given by 

e 2 (eV) 

X e ( m ) = SL (48) 

6.5 x 10“ 20 n e (m-3) 

This quantity is plotted in Fig. 36 as a function of plasma density with 
primary electron energy as a parameter. 

There are .some significant features of Fig. 36 that are noteworthy. 
Typical internal plasma conditions for both mercury and xenon cathodes are 
described by a plasma density of * 2 x 10“ cnf3 and a plasma potential (and 
hence primary electron energy) of ^9 v. Figure 36 shows these conditions 
result in an electron elastic scattering mean free path on the. order o' r 
6 mm. Because the energy exchange. mean free path for primary electrons 
(x pr ) goes as the sum of the inverses of A Jn and (Eq. 47), its numerical 
value is dominated by the shorter of these two mean free paths. For mercury 
x In is % 1 nim and this term dominates so A pr is ^ 1 mm. On the 

other hand for xenon, Aj n is typically much larger than- so that x pr is 
generally close in length to A 0 . Physically, this means that in the case 
of the mercury plasma the primary electrons give up most of their ener.gy 
in excitation and ionization reactions, while for xenon they tend to be 
thermal i zed by elastic scattering with Maxwellian electrons. This, suggests 
that for the xenon plasma much of the excitation and ionization is induced 
indirectly from the high energy tail of the Maxwellian electrons. Such. a 
situation is also consistent with the fact that the xenon discharge is 
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observed to operate. at plasma potentials close. to the first excitation level. 

It is also Interesting. to note that. the effects of primary electron . 
energy on the inelastic and elastic scattering mean free paths (x In and \ Q ). 
shown for xenon in Figs. 34. and 36 are opposite. Hence, while elastic mean 
free path dominates at low primary electron energies, the situation reverses 
at energies above ^ 12 eV. Above this primary electron energy the in- 
elastic reactions begin to dominate and the situation resembles that found 
in the mercury discharge where most of the primaries give up their energy 
in excitation and ionization reactions. In fact, such a shift probably 
does occur if the insert surface work function increases significantly, since 
the work function increase is accompanied by an increase in plasma potential 
(see Fig. 31 ). 

Comparison of Result s with Model P redic tions 
Data were collected for the hollow, cathode operating on xenon and. on 
argon over a range of discharge currents and flow rates. The results of 
these experiments will be presented here and compared with predicted results 
based on calculations made using the analytical model. 

The key points of the analytical model were described briefly in an 
earlier section of this report and are discussed in detail in Ref. 28. A 
summary of the equations used in the model are reproduced from Ref. 28 in 
Appendix B., If controllable operating parameters (discharge current and 
mass flow rate) .and the. cathode physical configuration (dimensions, surface 
work function and thermal characteristics are specified, the equations of 
the model, as. summarized in Appendix B, can predict cathode operating con^ 
ditions such as., insert temperature and plasma potential*. Calculations using 

the model were performed for both xenon and argon based on the following 

considerations. 
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1) Experimentally determined pressures, were used as Input to the. model 
rather than mass flow rate. In design applications where, measured pressures 
are not available, the empirical relationship given in Eq, B1 of .the model 
can be used to estimate the pressure from the flow rate and orifice diameter. 
The experimental results used in determining the parameters of the empirical 

pressure-flow rate relation (r . B1 ) for argon and xenon are. presented in 
the next section. 

2) The Surface work function $ s for the insert was not known. The 
value. used in the calculations was. chosen to give agreement between the 
calculated temperature and the average measured temperature at one arbitrarily, 
chosen experimental condition. The comparison between the calculated and 
measured average temperatures will, therefore, be a relative one valid mainly 

for checking the functional dependence on discharge current and pressure 
predicted by the model. 

3) The thermal power loss from the insert (Q th in Eq. B6) used for 
the calculations was based on the thermal analysis described in Appendix C 
of Ref. 30 for what is essentially the same insert-cathode tube config- 
uration used in the present experiments. The following equation represents 
an approximation to that analysis. 

Q th " (526 T^- 3.99 x 10 5 ) L g d in . (49) 

In this equation, T g is the average emission surface temperature in °K and 
the emission length L 0 and insert diameter d. p are both in meters. Equation 49 
was used in calculating the results presented here. 

4) The electron temperature was assumed to be 0.8 eV for the xenon 
calculations and 1.3 eV for argon. 

5) The ratio of emission length to mean free path L./a was chosen 

g pr 

to be one half in applying the model to argon and xenon. A ratio of unity 
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would indicate that the. active region of the cathode was indeed one primary 
energy_exchange mean free path long.. Such a criterion for the active region 
was, however, never expected .to be precise and the. ratio, of .one half used in 
the calculations was simply chosen to give reasonable agreement with the 
experimentally observed emission lengths. For mercury, a ratio of two was 
found to give good results. The fact, that a value of 0.5 gives good results 
for both argon and xenon may be related to the fact that the excitation cross- 
section data for both gases were drawn from the same source (Ref. 29). Other 
cross-section data for xenon, showed the excitation cross-section in- 
creasing almost twice as fast near threshold as that used here from Ref. 29. 

If the data of Ref. 31 were used, It would significantly alter the ratio 
of L e / A pr » pushing it closer to unity. 

6).... Because Eq. BIO cannot be solved explicitly for the surface 
temperature T s the equations in the model were solved in an iterative 
manner. Such a solution requires a determination of the total inelastic 
cross-section at each iteration. Since the calculations were made on a 
computer, this was handled by entering.the cross-section data from Fig. 33 
as arrays of data pairs and using an interpolation routine to determine 
intermediate values. 

The results of calculations using the model based on the assumptions 
discussed above are shown in Figs. 37 through 40 for argon and Figs, 41 
through 44 for xenon. The predictions of the model are plotted as the 
solid curves, while experimental data points corresponding to the predicted 
parameter are indicated by the circles. All of the experimental data were 
obtained using a cathode operating with a 0.76 mm diameter orifice. In each 
pair of plots, the upper plot shows the indicated parameter as a function 
of discharge current at a constant flow rate (287 mA for argon; 92 mA for 
xenon); while the lower plot Shows the parameter plotted as a function of 
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internal pressure at a constant discharge current of 2.3 A. The calculations 
for argon were based on an insert surface work function of 2.5 eV; while 
those for xenon were for a work function of 2.25 eV. 

The effects of discharge current, and internal pressure on emission 
surface length are shown in Fig. 37 for operation on argon. The experimental 
emission- lengths were determined from. the measured insert temperature 
profiles using the criterion that the emission extended upstream to a point 
where the temperature (in °K) had dropped to. 90% of the maximum temperature. 
It is estimated that a length based on this criterion accounts for 'v 90% of 
the total insert emission current. Figure 37 shows that the experimentally 
determined emission length increases moderately with increasing discharge 
current at a constant flow rate (upper plot) and is relatively independent 
of internal pressure (lower plot) over the range of pressures covered in 
the experiments. The calculated emission lengths based on a criterion of 
one half of the primary energy exchange mean free path are shown in Fig. 37 
by the . sol id curve. The sol.ld curves show discharge current and pressure 
trends similar to the. experimental data. 

At this point,. it is worth pointing out a general feature of the 
experimental results shown in. Fig. 37 which is also present in most of the 
other figures. In both the upper and lower plots of Fig. 37 all of the data 
fall along a single trend line except for one data point. The data along 
the trend line were all taken in sequence. For example, in the upper plot 
the discharge current was varied while holding the flow rate constant to 
obtain the data which follow the trend line. The odd point in that plot 
comes from the corresponding experiment used in obtaining the trend line 
data in the lower plot. That is the experiment in which the discharge 
current was held constant while the pressure was varied. The difference 
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between the odd date point and the trend line Is an indication of the change 
in operatlng.conditlons which took place between the_d 1 fferent phases of-the. 
experiment and Is believed to be due mostly to the changes in the Insert 
surface condition which were discussed earlier. 

Figure 38 shows the effect of discharge current and internal pressure 
on the average emission surface temperature for the cathode operating on 
argon. The average temperatures represented by the circles are an integrated 
average (over the emission length) determined from the measured Insert 
temperature profiles. Both the experimentally determined and the calculated 
surface temperature are seen 1n_F1g. 38 to increase with increases in dis- 
charge current and with Increases in internal pressure. As mentioned earlier, 
the calculated results were based on a surface work function of 2.5 eV chosen 
to give agreement at one arbitrary experimental condition (apparent in Fig. 38) 
The comparison between the calculated and experimental results is, therefore, 
on the basis of curve shape and slope. Based on this criterion the agreement 
between the calculated and experimental temperature results is reasonably 

good, particularly for the temperature as a function of internal pressure 
(lower plot). 

The Increase in emission temperature with Increasing pressure (lower 
plot. Fig. 38) seen with argon is just opposite to the mercury results which 
showed a decrease in temperature with Increased pressure . 28 In both cases 
(argon and mercury), the model predicts an. Increase in temperature of the 
magnitude and shape shown in Fig. 38. This predicted behavior Is linked to 
the decrease In emission area with increasing pressure (lower plot. Fig. 37 ). 

It is unclear why, in the case of mercury, the measured temperature did not 
also follow such a trend. . 

The effects of discharge current and internal pressure on the plasma 
density In the ion production region are shown In Fig. 39 for the cathode 
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operating on argoru. Again thejolld curve indicates the behavior predicted., 
by the. model. The experimental results_indicated by the circles were obtained 
from two different measurements. The open symbols and dashed line were based 
on the spherical probe measurement and are the maximum value of the plasma 
density along the cathode axis. The solid circles are results calculated 
from the measured current to the upstream surface of the orifice plate 
assuming that the current was due solely to the Bohm flux of ions to that 
surface. As such, the solid symbols represent what is really an average 
plasma density across the downstream boundary of the ion production region. 
Since the model assumes uniform properties within the ion production region, 
the plasma density calculated. with the model is really an average density 
and is more appropriately compared with the average measured values indicated 
by the solid symbols than with the maximum values-fopen symbols).. Such a 
comparison shows the predictions of the model to be quite good both in mag- 
nitude and functional dependence on the Independent variable. 

It is interesting to note here that the plasma density, both measured 
and predicted, are relatively constant as the pressure is increased at a 
constant discharge current (lower plot. Fig. 39). This is somewhat sur- 
prising since one would expect that to first order the plasma density would 
increase as the neutral density increases. The fact that this does not 
occur suggests competing effects such as reductions in excitation and ioni- 
zation cross-sections due to decreased electron energy. . This is, in fact, 
observed since both electron temperature (Fig. 32) and primary energy (see 

Plasma potentials, lower plot, Figs. 40 and 44) were observed to decrease 
with internal pressure. 

In Fig. 40, the effects of discharge current and internal pressure on 
anode, keeper, and Internal plasma potential are shown for argon operation. 
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The keeper and anode potential are highly dependent on electrode geometry 
and spacing and downstream conditions such as background pressure. . They are 

shown In Fig,. 40 -for two reasons. First, to give a general, indication of 

where the potential drop occurs in the cathode; andjsecond, .to allow a com- 
parison of similar results obtained with xenon using the same anode and 
keeper configuration. The parameter predicted by the model is the internal 
plasma potential and is indicated in Fig. 40 by the solid curves.. Comparison 
between those curves and the circular symbols suggests that the plasma »• 
tentials predicted by the model are reasonably close both in their magnitude 
and in their functional dependence on the indicated parameters. 

The model consistently predicts a plasma potential somewhat lower than 
the measured potential. This was also the case with mercury^ 8 and will be 
seen to be the case with xenon (Fig. 44) as well. The cause for the. con- 
sistently low predicted values is believed to be related to a term neglected 
in the energy balance on which, the plasma potential calculation is based. 

That term is the energy flux. from the ion production region associated with 
the resonant and metastable atoms which leave the region. With the present 
model there is no reasonable way of accounting for that term and, as the 
results indicate neglecting the term still results in good qualitative agree- 
ment between the. model and the experiment. 

Results. of the kind discussed above are shown in Figs. 41 to 44 for a 
similar cathode operating on xenon. It should.be noted that-the Insert used 
in the tests in which these results were obtained had a work function somewhat 
lower (2.25 eV) .than -that used. in the argon tests (2.5 eV). However, with 
one exception the results shown in Figs,. 41 to 44 are qualitatively very 
similar to those presented in Figs. 37 to 40 for argon and will not be dis- 
cussed in detail except to make a few brief points of comparison. 


The main difference to be pointed out Is the numerous inflection points 
and discontinuity in the slope. of the calculated results for' the constant 
flow rate case (upper plot. Figs. .41 to 44),. The inflections in the curves 
for the calculated results correspond to the inflections in the inelastic 
cross-section curve for xenon (Fig. 33) at low primary electron energies. 

The dashed portion of the curve following the discontinuity in slope (upper 
plot, Figs. 41 to 44) indicates the conditions where the model predicts 
primary electron. energies below the 8.3 eV excitation threshold of xenon. 
This is apparent if one examines the predicted plasma potential curve in 
the upper plot of Fig. 44. At discharge currents above 'v 3.3 A, the model 
predicts an internal plasma potential less than 8.3 eV. The discontinuity 
in the slope of the plasma potential curve at a discharge current of 3.3 A 
is, therefore, just a reflection of the discontinuity in the slope of the 
cross-section curve at the threshold energy of 8.3 eV. Physically, the 
solutions, at .discharge currents above 3.3 A represent the situation where 
most of the excitation and -ionization must be due to the energetic tail of 
the Maxwellian electron distribution since the primaries have insufficient 
energy to excite or ionize ground state atoms. (Excitation and ionization 
of excited neutrals would still be possible.) This is the extreme case of 
the situation discussed earlier where it was noted that for xenon at lower 
primary energies the elastic scattering term for the primaries dominates 
over the inelastic one in determining the primary energy exchange mean free 
path. In such a situation most of the primary energy would be deposited in 
the Maxwellian electron population, and it is that population which must 
sustain the discharge. 

It is worth making a couple of additional comparisons between the 
xenon and argon results. Both the predicted and the measured values of 
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emission length for xenon (Fig. 41) were substantially shorten than for 
argon (Fig. 37)... The average plasnm_density. is greater for xenon. than for 
argon (lower plot* Figs, 39 and 43.). Finally, the anode potentials with 
xenon (Fig., 44) are 8 - 10 volts less than with argon (Fig. 40), 

Pressure-Flow Rate Correlations 

In the experiments just described, the. internal pressure was determined 
by direct measurement using the manometer shown in Fig. 28. However, 
normally the internal cathode pressure is unknown and mass flow rate is the 
parameter controlled during cathode operation. For most applications of 
the model, it is, therefore,. .desirable to have some means of relating in- 
ternal cathode pressure to mass flow rate. Predicting this pressure from 
cathode dimensions,. mass flow rate, and.operating conditions is, however, .a 
complex problem made particularly difficult by the following effects. ..For 
normal cathode dimensions and flow rates, the cathode operates in the transi- 
tion regime between free molecular and continuum flow. In addition, because 
the gas is ionized, the flow is subject to energy input from plasma processes 
and is affected. by the forces due to the electric fields set up within the 
plasma. Because of the complexity of the problem, the model developed for 

cathodes operating on mercury used an empirical formula for predicting 

28 

cathode internal pressure. 1 The same approach will be taken here for 
argon and xenon. 

In order to determine an empirical relation between th e p ressure and 
the flow rate, an experiment was run in which the stagnation pressure in 
the plenum chamber immediately upstream .of the cathode was measured with 
the manometer for a wide range of emission currents (0-4. 3A), orifice diam- 
eters (0.25 to 1.0 mm), and mass flow rates (50 to 1300 mA). The throttle 
valve shown in Fig. 28 was, of course, closed during these tests. The results 
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are shown In Figs. 45 and 46 respectively for argon. and xenon where the 
parameter P 0 /(m/d 2 ) is plotted as a function. of total discharge current. 

Here P Is the total stagnation pressure in Torr, m is the mass flow rate 
In mi 11 i amps equivalent, and d Q is the orifice diameter in millimeters. 

The circles in Figs. 45 and 46 represent the average value for all of the 
data collected at that discharge current. A straight line can be used to 
fit the average values reasonably well . .The .aquations shown in Figs. 45 
end 46 represent the least squares fit of the data points to a straight 

line. 

The bar on the symbols in Figs. 45 and 46 represents the range of the. 
data for the various orifice sizes and flow rates. The full range is shown 
by the bars, rather than the standard deviation because the spread in the 
data is not believed to be due so much to random experimental error but to 
the fact that the form of the correlation factor P Q /(m/d 2 ) does not com- 
pletely account for the influence of the flow rate and the orifice diameter. 
It is believed that if the correlation factor fully accounted for the param- 
eters of the problem, the range indicated by the bars would be much smaller. 
This form of the correlation factor was chosen, however, because it was used 
earlier, in. the mercury work 28 and it has a theoretical basis. Both the 
theory of free molecular flow and the theory of continuum, choked flow pre- 
dict that the pressure-mass flow relation takes the form 

P 0 /(m/d 2 ) = C/T7M (50) 

where T is the stagnation temperature at the orifice and M is the atomic... 
mass of the gas. The parameter C is a constant of proportionality which is 
different For the two cases (free molecular and continuum). The gas tern 
perature for the zero discharge. current case was maintained during the tests 
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Fig. 46. Internal Pressure/ Flow Rate Correlation for Xenon 
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at. 1 '' 1Q00°C using the. cathode heater. (This is probably also close, to the 
neutral particle temperature when the discharge is on).. The values of 
c'/ffi based on this temperature are plotted on the vertical axis for the 
free molecular and continuum flow cases. For the no discharge case, the 
hollow cathode flow for both argon and for xenon falls in the transition 
region between the two types of flow, though it is much-closer to being 
continuum flow than free molecular flow. The parameter c'/TM is used here 
to distinguish it from.. the term on the right hand side of Eq. 50. The . 
parameter c'/TM accounts for the fact that the units of m used for the 
correlation factor in Figs, 45 and 46 are milliamps equivalent. That is 
m as ..used in the figures really represents a. .charged particle flow. rate 
rather than a mass flow rate, hence,. ..the dependence /fM rather than /f/H,. 

28 

For comparison, the pressure flow correlation determined for mercury 
is shown in Figs., 45 and 46_hy a dashed line. The center dashed line shown 
in Fig. 45 was obtained by multiplying. the mercury result-by 
y/M ( argon ) /M (me rcu ryj. Predicting argon behavior from the mercury result in 
this manner is seen to yield good agreement in the no discharge case. How- 
ever, when the discharge is present such an extrapolation of the mercury 
data predicts a pressure higher than the measured argon pressure. Similar 
results are shown in Fig. 46 for xenon. While it was initially anti ci <- .eu 
that the pressure flow. results should scalawith atomic mass in the manner 
suggested by C /fM, this is obviously not the case when a discharge is 
present. 

It . is suggested that the equations indicated in Figs. 45 and 46 be 
used to calculate the internal pressure when applying the model to cathodes 
using argon or xenon and the pressure is not known. This is indicated by 
Eq. Bl which appears in the summary of the model presented in Appendix B. 
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Conclusions 

The .results_of the_ tests discussed above suggest that the physical 
processes important in the operation of argon and xenon hollow cathodes are 
similar to those for the mercury cathode. Namely, that field-enhanced 
thermionic emission is the predominant surface emission mechanism with a 
significant fraction. of the total current coming from the Bohm flux of ions 
to the cathode surfaces. Surface emission occurs mainly from a region near 

the downstream end of the insert and extends over an axial distance on the 

order of the primary electron energy exchange mean free path. The one 
notable difference between the processes found in the argon and xenon ca- 
thodes and those in the mercury cathode is that much of the excitation and 
ionization is carried out by the high energy tail of the Maxwellian electrons 
for the inert gases. This effect is most significant at low insert surface 
work function conditions where the plasma potential is sufficiently low so that 
the primary energy drops below the excitation threshold energy. In this 
situation, the inelastic contribution to the primary electron. energy exchange 

mean free path becomes very small and most of the primary electrons are 

thermalized before having an inelastic collision. However, this difference 
in the collisional processes in the cathode_does net invalidate the analytical 
model. In fact, when applied to hollow cathodes operating on argon and 
xenon, the model was found to predict operating parameters which agreed 
well with the experimental results both_in magnitude and in their functional 
dependence on operating conditions. This agreement was predicated on the 
assumption that the ratio of emission length to primary electron energy 
exchange mean free path was one half. This is a reasonable value given the 
general assumptions and approximations made in the model. The fact, that 
predictions of the model agree well with experimental results for three 
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propellants which have significantly different excitation cross-sections 
and atomic masses is a strong indication that the assumptions of the model 
are valid and represent a good description of the basic physical processes 
underlying the operation of the orificed hollow cathode. 

For a given insert surface work function the emission surface tempera- 
tures are similar for cathodes operating jsn all three propel! ants--argon, 
xenon* and mercury. However, for operation with argon and xenon low sur- 
face work functions were difficult to obtain and they tended to increase 
significantly during the period of several hours required to conduct the 
experiments. It is believed that the degradation of the. insert is due to 
prevailing plasma conditions which enhance sputter damage of the surface 
when operating with argon and xenon. As the surface work function increases 
the surface temperature and plasma potential both increase causing increas- 
ingly rapid deterioration of the surface. It is believed that this degrada- 
tion process is probably significantly worse with the single layer tantalum 
foil insert used in these experiments than it would be with the impregnated 
inserts used in actual thruster cathodes. It is quite possible that the 
surface work function of impregnated inserts would appear to be relatively 
stable while increased sputtering of the surface would increase the depletion 
rate of the impregnant and, thereby, reduce the lifetime of the cathode. 

This is an area that should probably be investigated more closely if cathodes 
with impregnated inserts are to be operated with argon and xenon. 

The present experiments suggest that rolled foil inserts coated with 
a double carbonate mixture may not exhibit adequate lifetimes when operated 
with argon and xenon. An investigation into alternative insert materials 
and/Or coatings or the use of a multilayer foil insert without any coating 
would also be indicated as potentially fruitful areas for improving 


performance and reliability of cathodes operating with argon and xenon. 
Aston and Kaufman have both had some success with uncoated -tantalum 
foil inserts. Such a cathode, however, would require careful design to 
minimize thermal losses from the insert and special provisions to insure 
reliable starting. In addition, discharge voltages in such a cathode are 
expected to be significantly higher than. in a cathode with a coated or 
impregnated insert. 
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DISCHARGE CHAMBER WORK 


' A g - Grid area (m 2 ) 

; c o - 4 Ve /(ev o A g »o> 

I e- - Electronic charge (1.6 x 10" 19 coul . ) 

f B - Extracted ion fraction 

■< f c - Fraction of ion current produced that goes to cathode 

r u potential surfaces. 

; 0^ - Ion current to anode potential surfaces (A) 

] Jg - Beam current (A) 

,-x 

% J c -Ion current to cathode potential surfaces (A) 

i J D - Discharge current (A) 

Jg - Cathode emission current (A) 

J* L. 

J. - Rate of production of j tn excited state expressed as 
J a current (A) 

0 L - Current of primary electrons lost to the anode (A) 

J^ - Current of Maxwellian electrons to the anode 

Jp - Ion production current (A) 
k - Boltzmannls constant (1.38 x 10‘ 23 J/°K) 

_ & e - Primary electron containment length (m) 

m^ r - Argon atomic mass (kg) 

\ r - Krypton atomic mass (kg) 
m - Propellant mass flow rate (Aeq.) 
n fi - Plasma density (nr 3 ) 
n Q - Neutral density (m" 3 ) 
n Q - Neutral loss rate (Aeq.) 
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- Probability of a primary electron/neutral atom collision 
prior to electron collection at the anode 

- Neutral atom temperature (°k) 

- Wall temperature (°k) 

- Excitation energy of j th excited state (eV) 

- Ionization energy (eV) 

- Potential (relative to cath. pot.) from which electrons are 
accelerated to become primaries (V) 

- Discharge voltage (V) 

- Electron velocity (m/s) 

- Neutral atom velocity (m/s) 

- Plasma volume (m 3 ) 

- Average beam ion energy cost (eV) 

- Average energy of Maxwellian electrons leaving the plasma at 
the anode (eV) 

- Average plasma ion energy cost (eV) 

- Average plasma ion energy cost neglecting primary electron 
energy lost .to the anode (eV) . 

- Average plasma ion energy cost considering ionization and 
excitation processes only (eV) 

- Propellant utilization 

- Total inelastic collision cross section for primary electron - 
neutral atom collisions. (m 2 ) 

- Ionization collision cross section (m 2 ) 

- Excitation collision cross section of j th state (m ) 

- Effective grid transparency to neutrals 

- Accelerator grid transparency 

- Screen grid transparency 
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APPENDIX B 

EQUATION SUMMARY FOR HOLLOW CATHODE MODEL 
Dan Siegfried 

The collection of equations which comprise the analytical model des- 
cribing hollow cathode operation are shown in Table BI. Originally 

developed for a cathode operating on mercury, these equations appear to hold 
equally well for operation on argon and xenon, and are the ones used to . 
calculate the results presented. The first column of Table BI. indicates 
the physical basis for the equation shown in column three and the second 
column notes important assumptions. Table BI I. gives the empirically de- 
termined constants (C lv and C 2 ) used in the pressure-flow rate correlation 
given by Eq. BI . and also indicates recommended values for the constant 
C 3 and the electron temperature (T g ). Table Bill, lists the symbols and 
units used for the parameters in the equations in Table BI. It should be 
noted that all of the equations in Table BI. are in MKS units EXCEPT for the 
parameters in the empirical expression for the pressure (Eq. BI). 

Given the mass flow rate m, the orifice diameter d 0 , the total dis- 
charge current Ip , the surface work function <f> $ , the insert thermal power 
loss Q th , and the physical constants (e, k, m.. , e Q , e^, a Q ), the equations 
in Table BI. can be solved for all of the other parameters except the 
electron temperature. Suggested values for the electron temperature are 
indicated in Table BI I . Because Eq. BIO cannot be solved explicitly for 
the surface temperature T $ , the solution of these equations is necessarily 
iterative. Such a solution requires a determination of the total inelastic 
cross-section at each iteration. If the calculations are made on a com- 
puter, this can be handled by entering the cross-section data as an array 
of data pairs and using an interpolation routine to determine intermediate 
values. - - - 


Summary of Equations Used in Model 
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Table 6,111 

NQMANCLATURE FOR CATHODE WORK 

a Q - Theoretical constant ( 1 .2x1 0 6 A/m 2 °K 2 ) 

A c - Area of end boundary of ion production region. (in 2 ) 

A e - Insert emission area (m 2 ) 

A s - Total -surface area of ion production region (m 2 ) 

Cj - Empirical constant for pressure-flow rate correlation 
(SeeTableB.il.) 

C 2 - Empirical constant for pressure-flow rate correlation 
(See Table B.II.) 

C 3 - Empirical constant relating emission length and mean free path 
(See Table B.II.) 

d Q - Orifice diameter (mm) 

e - Electronic charge {Coulombs) 

E - Electric field at insert surface (V/m) 

I D - Total discharge current (A) 

I g - Insert electron emission current (A) 

1^ - Total ion current to cathode surfaces (A) 

J\j - Bohm current density (A/m 2 ) 

J th - Field-enhanced thermionic emission current density (A/m 2 ) 
k - Boltzman's constant (1 .38xlO" 23 J/°K) 

Lg - Insert emission length (m) 
m - Propellant mass flow rate (mA equivalent) 
m - Ionic mass (kg/ion) 

n 0 - Electron density in ion production region (m -3 ) 

n Q " Total neutral acorn density in ion production region (nf 3 ) 
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Table B.III (continued) 

P - Internal cathode pressure 
Q th - Insert thermal power loss (W) 

Tp “ Electron temperature in ion production region (°K) 

T $ - Insert surface temperature (°K) 

Vp ~ Plasma potential in ion production region (V) 

e.. - Ionization potential 

e Q - Permittivitity of free space 

e pr “ Primary electron energy, equivalent to plasma potential V (eV) 

r 

x pr " Primar y electron, energy exchange mean free path (m) 

<t> e - Average effective work function (V) 

<|> s - Average surface work function (V) 


